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Abstract 

Two long-standing problems with the post-Newtonian approximation for iso¬ 
lated slowly-moving systems in general relativity are : (i) the appearance at 
high post-Newtonian orders of divergent Poisson integrals, casting a doubt on 
the soundness of the post-Newtonian series; (ii) the domain of validity of the 
approximation which is limited to the near-zone of the source, and prevents 
one, a priori, from incorporating the condition of no-incoming radiation, to 
be imposed at past null infinity. In this article, we resolve the problem (i) 
by iterating the post-Newtonian hierarchy of equations by means of a new 
(Poisson-type) integral operator that is free of divergencies, and the prob¬ 
lem (ii) by matching the post-Newtonian near-zone field to the exterior field 
of the source, known from previous work as a multipolar-post-Minkowskian 
expansion satisfying the relevant boundary conditions at infinity. As a re¬ 
sult, we obtain an algorithm for iterating the post-Newtonian series up to 
any order, and we determine the terms, present in the post-Newtonian field, 
that are associated with the gravitational-radiation reaction onto an isolated 

slowly-moving matter system. 
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I. INTRODUCTION 


A. Problems with the post-Newtonian expansion 


The post-Newtonian approximation, or expansion when the speed of light c —> +cx), has 
been formalized in the early days of general relativity by Einstein |l[], Droste 0, and DeSitter 
0. Since then, it has provided us with our best insights about the problems of motion and 
gravitational radiation, two of general relativity’s most important issues. Concerning the 
problem of motion, we quote the dynamics of N separated bodies at the first post-Newtonian 
(IPN, or 1/c^) order : works of Einstein, Infeld and Hoffmann ^ and other authors |^|-0, 
and the dynamics of extended fluid systems up to the 2.5PN level of gravitational radiation 
reaction : works of Chandrasekhar and collaborators [^|-0 and followers In the 

case of two compact objects, we know the 2.5PN equations of motion of the binary pulsar 
, and the 3PN equations of motion of inspiralling compact binaries . The spe¬ 

cific contribution of the gravitational-radiation reaction has been obtained up to the 1.5PN 
relative order by the method of matched asymptotic expansions for extended fluids 


and by means of balance equations for compact binary systems Concerning the 

problem of gravitational radiation, the work has focused on the expressions of the multipole 
moments of general fluid systems BH. and on the gravitational-wave flux emitted by 
inspiralling compact binaries, including the specific effects of wave tails, up to the 3.5PN 
order 

The “standard” post-Newtonian approximation, at the basis of most of the body of 
work quoted previously, is known to be plagued with some apparently inherent difficulties, 
which crop up at some high post-Newtonian order like 3PN. Up to the 2.5PN order the 
approximation can be worked out without problems, and at the 3PN order the problems 


can be solved specifically for each case at hands (see for instance Ref. [^|). However, it 
must be admitted that these difficulties, even appearing at higher approximations, cast 
doubt on the actual soundness, on the theoretical point of view, of the post-Newtonian 
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expansion. What is maybe worse, they pose the practical question of the reliability of 
this approximation when comparing the theory’s predictions with very precise experimental 
resnlts. It is therefore highly desirable to assess the natnre of these difhcnlties - are they 
pnrely technical or linked with some fnndamental drawback of the approximation scheme? 
- and eventnally to resolve them. This is especially important in view of the fact that 
inspiralling compact binaries, when they are detected and analyzed by gravitational-wave 
experiments, will necessitate a prior theoretical knowledge of the gravitational-wave signal at 
some very high post-Newtonian order i4l|-^ . In this article let us distinguish (and resolve) 
the two basic problems faced by the post-Newtonian expansion. 

The first problem is that in higher approximations some divergent Poisson-type integrals 
appear. Recall that the post-Newtonian expansion replaces the resolution of an hyperbolic- 
like d’Alembertian equation by a perturbatively equivalent hierarchy of elliptic-like Poisson 
equations. Rapidly it is found during the post-Newtonian iteration that the right-hand-side 
of the Poisson equations acquires a non-compact support (it is distributed over all space), 
and that the standard Poisson integral diverges because of the bound of the integral at spatial 
inhnity, i.e. r = |x| —-|-cx), with t = const. For instance some of the potentials occuring 
at the 2PN order in Chandrasekhar’s work [Q are divergent, so the corresponding metric is 
formally inhnite Q . In fact, Kerlick showed that the post-Newtonian computation 

d la Chandrasekhar M. following the iteration scheme of Anderson and DeCanio [j^ 


can be made well-dehned up to the 2.5PN order, by keeping some derivatives inside some 
crucial integrals to make them hnite However, the latter remedy does not solve the 

problem at the next 3PN order, which has been found to involve some inexorably divergent 
Poisson integrals [0,^. 

These divergencies come from the fact that the post-Newtonian expansion is actually a 


^Nevertheless, these divergencies were not a problem when considering the equations of motion 
because the gradients of these potentials, which parametrize the equations, were finite. 
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singular perturbation, in the sense that the coefficients of the successive powers of 1/c are 
not uniformly valid in space, since they typically blow up at spatial inhnity like some pos¬ 
itive powers of r. For instance, Rendall has shown that the post-Newtonian expansion 
cannot be “asymptotically flat” starting at the 2PN or 3PN level, depending on the adopted 
coordinate system. The result is that the Poisson integrals are in general badly-behaving 
at inhnity. Physically this can be understood by the fact that the post-Newtonian approx¬ 
imation is valid only in the near zone of the source (see below) while the Poisson integral 
extends over the whole three-dimensional space, including the regions far from the source 
where the approximation breaks down. Therefore, trying to solve the post-Newtonian equa¬ 
tions by means of the standard Poisson integral does not a priori make sense. This does not 
mean that there are no solution to the problem, but simply that the Poisson integral does 
not constitute the correct solution of the Poisson equation in the context of post-Newtonian 
expansions. So the difficulty is purely of a technical nature, and will be solved once we 
succeed in hnding the appropriate solution to the Poisson equation Q. A solution to the 
problem of divergencies has been proposed by Futamase and Schutz and Futamase 


Their approach is alternative to the one we shall follow below. It is based on an initial- 
value formalism, which avoids the appearance of divergencies because of the hniteness of the 
integration region. 

The second problem has to do with the near-zone limitation of the approximation. Indeed 


^The problem is somewhat similar to what happens in Newtonian cosmology. Here we have to 
solve the Poisson equation AC/ = —AnGp, where the density p of the cosmological fluid is constant 
all over space : p = p{t). Clearly the Poisson integral of a constant density does not make sense, as 
it diverges at the bound at infinity like the integral f rdr. This nonsensical result has occasionally 
been referred to as the “paradox of Seeliger”. However the problem is solved once we realize that 
the Poisson integral does not constitute the appropriate solution of the Poisson equation in the 
context of Newtonian cosmology. A well-defined solution is simply given by U = — |7rGpr^. 
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the post-Newtonian expansion assumes that all retardations r/c are small, so it can be viewed 
as a formal near-zone expansion when r —> 0, which is valid only in the region surrounding 
the source that is of small extent with respect to the typical wavelength of the emitted 
radiation : r <C A (if we locate the origin of the coordinates r = 0 inside the source). 
Therefore, the fact that the coefficients of the post-Newtonian expansion blow up at spatial 
infinity, when r —>■ -|-cx), has nothing to do with the actual behaviour of the field at infinity. 
The serious consequence is that it is not possible, a priori, to implement within the post- 
Newtonian iteration the physical information that the matter system is isolated from the 
rest of the universe. Most importantly, the no-incoming radiation condition, imposed at past 
null infinity, cannot be taken into account, a priori, into the scheme. In a sense the post- 
Newtonian approximation is not “self-supporting”, because it necessitates some information 
taken from outside its own domain of validity. 

To the lowest post-Newtonian orders one can circumvent this difficulty by considering 
retarded integrals that are formally expanded when c —^ -1-cxd as series of “instantaneous” 
Poisson-like integrals |^. This procedure works well up to the 2.5PN level and has been 
shown to correctly fix the dominant radiation reaction term at the 2.5PN order |T^,IT3 


Unfortunately such a procedure assumes fundamentally that the gravitational field, after 
expansion of all retardations r/c —> 0, depends on the state of the source at a single time, 
in keeping with the instantaneous character of the Newtonian interaction. However, we 
know that from the 4PN order the post-Newtonian field (as well as the source’s dynamics) 
ceases to be given by a functional of the source parameters at a single time, because of the 
imprint of gravitational-wave tails in the near zone field, in the form of some modification, 
at the 1.5PN relative order, of the radiation reaction force |3T|-|33|. Therefore, the formal 
post-Newtonian expansion of retarded Green functions is no longer valid starting at the 
4PN order. We face here a true difficulty, which is fundamentally linked to the nature of 
the post-Newtonian approximation. 

The aim of the present article is to resolve the two latter problems. We shall prove 
that the post-Newtonian expansion can be indefinitely reiterated, while incorporating the 
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correct boundary conditions satisfied by the wave field at infinity. In particular, we shall 
get new insights about the problem of gravitational-radiation reaction inside an isolated 
(post-Newtonian) system. To cure the problem of divergencies we introduce, at any post- 
Newtonian order, a generalized solution of the Poisson equation with non-compact support 
source, in the form of an appropriate finite part of the usual Poisson integral : namely we 
regularize the bound at infinity of the Poisson integral by means of a process of analytic 
continuation, analogous to the one already used to regularize the retarded integrals in Refs. 
i3^,|39|,^. Our generalized solution constitutes a particular (well-defined) solution of the 
problem; the most general solution is the sum of that particular solution and the most general 
solution of the corresponding homogeneous equation, i.e. the source-free Laplace equation. 
The homogeneous solution should be regular all over the matter system (we are considering 
smooth matter distributions), and we find, after summing up the post-Newtonian series, 
that it can be thoroughly written with the help of some tensorial functions of time Af^{t), 
where L = R-'-p denotes a multi-index with I indices At this stage, considering 

the post-Newtonian iteration scheme alone, we cannot do more and therefore we leave the 
functions Aff {t) unspecified. We refer to them as some “radiation-reaction” functions. 

The solution of the problem of the near-zone limitation of the post-Newtonian expansion 
resides in the matching of the near-zone field to the exterior field, a solution of the vacuum 


equations outside the source which has been developed in previous works p6|,|3^ using some 
post-Minkowskian and multipolar expansions. In the case of post-Newtonian sources, the 
near zone, i.e. r -C A, covers entirely the source, because the source’s radius itself is such that 
a -C A. Thus the near zone overlaps with the exterior zone where the multipole expansion is 
valid. Matching together the post-Newtonian and multipolar-post-Minkowskian solutions in 
this overlapping region is an application of the method of matched asymptotic expansions. 


and has frequently been applied in the present context, both for radiation-reaction [p9|-|33 
and wave-generation [p7|-|i0[] problems. 

The exterior multipolar-post-Minkowskian field originally obtained in Ref. depends 
on some “multipole-moment” functions, say Xff' {t) [whose components are associated with 
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some source multipole moments e.g. //.(t), Jhif), ■ ■ ■ ], which must be left unspecihed as long 
as we consider only the external vacuum solution. In the work , we have shown that the 
multipole moments {t) are entirely determined, up to any post-Newtonian order, from 
the requirement of matching to a post-Newtonian solution. In the present paper, we shall 
further show that the radiation-reaction functions parametrizing the post-Newtonian 

solution, are also uniquely hxed, up to any post-Newtonian order, by the matching. In 
particular, we shall hnd that the latter functions include correctly the contribution of wave 
tails, arising at the 4PN order, as determined in Refs. |l3l|-0^. We shall also recover by a 
different method the result of Ref. m concerning the multipole moments X^{t). 

A comment is in order regarding the possibility of determining the near-zone held by 
matched asymptotic expansions up to any post-Newtonian order. Indeed the method pre¬ 
supposes the existence of the exterior near-zone for which a < r X. Now if a given physical 
system, whose dynamics is described by Newton’s theory, emits gravitational radiation at 
some Newtonian fundamental wave length An, we expect that when taking into account 
the post-Newtonian corrections up to the post-Newtonian order n, it will have a radiation 
spectrum composed of harmonics between ~ 2AN/'n and ~ 2 An. Indeed this is the case of 
the radiation from a binary system moving on a circular orbit, for which we have ;^An < 
AnPN < 2 An. Therefore, if n is large enough, say n > 2AN/a, we expect that there will be 
some part of the radiation whose frequency is too high for the exterior near zone to exist. 
What we want to say is that the formulas we shall obtain for the post-Newtonian held of 
a source “up to any order” are indeed physically valid, strictly speaking, only up to some 
hnite post-Newtonian order ~ 2X^/a, where a is the size of the source; but that, if we 
consider a source which is less relativistic, for instance which is obtained by “slowing down” 
our source so that its Newtonian wave length gets twice its original value (say), the same 
post-Newtonian formulas can then be used for the new source up to approximately twice 
the previous post-Newtonian order. 

The plan of this paper is as follows. In Section we recall the construction in Ref. 

3^ of the multipole expansion of the external held, and we obtain thanks to a result of 
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Ref. the near-zone expansion of that external held ready for subsequent matching. In 
Section EH we implement the post-Newtonian iteration of the inner held inside the matter 
source, and we hnd the far-zone (multipolar) expansion of that post-Newtonian solution, 
also ready for matching. In Section 0 we show that the matching works up to any post- 
Newtonian order, and permits the determination of all the unknowns, in both the external 
and inner helds. Finally in Section 0 we check that our post-Newtonian solution satishes the 
harmonic-coordinate condition as a consequence of the equations of motion of the source. 
The technical proofs are relegated to Appendices 0, H and p. 


B. Notation for the Einstein field equations 

For the problem at hands let us introduce an asymptotically Minkowskian coordinate 
system for which the basic gravitational-wave amplitude, = y/—g is diver¬ 
genceless, i.e. satishes the de Bonder or harmonic gauge condition = 0. Here, g^'^ 

denotes the contravariant metric (satisfying g^^gpu = , g is the determinant of the covari¬ 

ant metric, g = det{gpt^), and represents an auxiliary Minkowskian metric with signature 
+2. With these dehnitions the Einstein held equations can be recast into the d’Alembertian 
equation 


= 


IGttG 


.fliy 


( 1 . 1 ) 


where □ = -|- A is the (hat-spacetime) d’Alembertian operator. The 


source term, can rightly be interpreted as the “ehective” stress-energy pseudo-tensor 
of the matter and gravitational helds in harmonic coordinates. It is conserved in the usual 
sense, and that is equivalent to the condition of harmonic coordinates : 




The pseudo-tensor is made of the contribution of the matter helds, described by a stress- 
energy tensor and the one due to the gravitational held, given by the gravitational source 
term thus, 


t"" = IjIT"- + . 


(1.3) 


The conservation property (O) is equivalent to the conservation, in the covariant sense, of 
the matter tensor : = 0. The exact expression of taking into account all the 

non-linearities of the Einstein held equations, reads 


am- = _ + -g^’^gpM^^drh 


crA 


+ ^{2g^Pg'^- - g^'^gni^gxrg.. - gregx.)dph^^. 


(1.4) 


It is clear from this expression that is made of terms which are at least quadratic in the 
gravitational-held strength and its hrst and second space-time derivatives. 

In this article, we look for the solutions of the held equations (|1.1| )- (|1.4| ) under the fol¬ 
lowing hypotheses. First, we assume that the matter tensor has a spatially compact 
support, i.e. can be enclosed into some time-like world tube, say r < a, where r = |x| is the 
harmonic-coordinate radial distance. Second, we assume that the matter distribution inside 
the source is smooth : i.e. T^'^(x, t) G We have in mind a smooth hydrodynamical 

“huid” system, without any singularities nor shocks (a priori), that is described by some 
Eulerian-type equations including high relativistic (post-Newtonian) corrections. In partic¬ 
ular, we exclude from the start any sources containing black holes. Notice, however, that 
it makes sense to apply the formulas derived a priori only for smooth matter distributions 
to systems containing compact objects (including black holes), described by some sort of 
point-particle singularities; see e.g. Refs. |4I|-^. Finally, in order to select the physically 
sensible solution of the held equations, we choose some boundary conditions at inhnity cor¬ 
responding to the famous no-incoming radiation condition. In this paper, we shall rely on 
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a specific construction of the metric outside the domain of the source (r > a), that was 
achieved in Ref. under the assumption that the gravitational held has been independent 
of time (stationary) in some remote past, in the sense that t < —T ^ t)] = 0. 

This condition is a mean to impose, by brute force, the no-incoming radiation condition R 


II. EXTERIOR FIELD 


A. Multipolar expansion of the non-linear vacuum field 


In this section we review some material from Ref. concerning the construction of 
vacuum metrics by means of mixed multipolar and post-Minkowskian (MPM) expansions. 
The so-called MPM metrics aim at describing the gravitational held in the region exterior 
to a general isolated system. In fact they are mathematically dehned in the open domain 
X R, i.e. R^ deprived from the spatial origin r = |x| = 0, but of course do not agree 
physically with the real solution when 0 < r < a, since they are vacuum solutions. For 
our present purpose the point is that the most general physically admissible solution of 


the vacuum held equations has been obtained in Ref. |^| by a specihc construction of the 
post-Minkowskian solution, say 


+ 00 




( 2 , 1 ) 


m=l 


whose coefficients are in the form of multipolar series, or equivalently decompositions in 
symmetric-trace-free (STF) products of unit vectors hi, that are equivalent to the usual 
decomposition in spherical harmonics : 


^However the condition of stationarity in the past, though much weaker than the actual no¬ 
incoming radiation condition, does not seem to entail any physical restriction on the applicability 
of the formalism, even in the case of sources which have always been radiating. 
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( 2 . 2 ) 


+ O0 


Vm > 1 , (x, t) = J2^L {0, 0) (r, t) . 


1=0 


The are certain functions of the radial coordinate r and of time t. Inserting the 

MPM expansion (|2.1|) - (|2.2|) into the vacuum held equations (|1.1|) - (|1.2|) we obtain, at any 
post-Minkowskian order m, 




(2.3) 

(2.4) 


where denotes the m-th post-Minkowskian piece of the gravitational source term dehned 
by Eq. (|1.4|) , i.e. in which we have inserted the previous post-Minkowskian iterations up 
to the previous order m — 1 [with the convention that = 0]. Because Eq. ( |1.4|) is at 
least quadratic in non-linearities, it is clear that only the preceding iterations, < m — 1, are 
necessary at any post-Minkowskian order m. 


Now the solution that was obtained in Ref. 36 has two main characteristics. The hrst 


one is related to its particular near-zone structure, which will play a fundamental role in the 
present article. Namely, it was proved that each one of the multipolar-post-Minkowskian 
coefficients in Eq. (|2.1|) - that we recall are only dehned when r > 0 -, admits a singular 
near-zone expansion, i.e. when r —0, owning the following structure : 


ViV e N, 





(2.5) 


where the multipolar order I G M, where the powers of r are such that a G Z with Omin < a < 
N (with Omin a negative integer), and where the powers of In r are p G N with p < m — 1 . The 
maximal divergence when r —> 0 occurs for amin, which depends on the post-Minkowskian 
order m, and satishes amin(^^) —^ —oo when m —>■ -|-cx). Similarly the maximal power 
of the logarithms, Pmax(^^) = m — 1, tends to inhnity with m. The functions ^^{t) 

are smooth functions of time, G C'°°(M), which are to be computed by means of 
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the algorithm proposed in Ref. [0, and appear as complicated non-linear functionals of 
some more elementary functions parametrizing the linearized (m = 1) approximation. The 
remainder term in Eq. (|2.5|) is such that 


t) = 0{r^) when r —> 0 and t = const 


( 2 . 6 ) 


The Landau (9-symbol takes its usual meaning. This remainder admits also some specihc 


differentiability properties (refer to for the details). The gravitational source term 
admits exactly the same near-zone structure as in Eq. (|2.5|) with the exception that 
Pmax = m — 2 in this case (that is, the maximal power of the logarithms increases by one 
unit when going from the source to the solution). 

The second important characteristic of the MPM solution concerns the constructive 
formula which dehnes it. We hnd that each one of the post-Minkowskian coefficients is 
explicitly constructed by means of the following : 


= FP DrI 

(m) B=0 


(m) 


H-(X) 


(t- " 


1=0 


(2.7) 

r ' 7- V Q/ 

The hrst term involves a special type of generalized inverse d’Alembertian operator, built 
on the standard retarded integral. 


which 

factor 


□ 


-1 

Ret 










( 2 . 8 ) 


|x-y| 

extends over the whole three-dimensional space, but inside which a regularization 
has been “artihcially” introduced, namely 


r^ = 



(2.9) 


where B denotes a complex number, i? G C, and tq represents an arbitrary constant length 
scale. The indication FP stands for the finite part at R = 0, and means that one should 
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first compute the Laurent expansion when B ^ 0 oi (the analytic continuation of) the 
i?-dependent integral (p.8|) , and, second, pick up the hnite part at 5 = 0 in that expansion, 
i.e. the coefficient of the zero-th power of B. The main property of this generalized retarded 
operator, that we shall from now on abbreviate as 


□ 


Ret 


FP □ 

B=0 


-1 

Ret 



( 2 . 10 ) 


is that, for source terms admitting a near-zone structure of the type (|2.5|) , 


□ 


□ “1 


- 

- 


( 2 . 11 ) 


Because the second term in Eq. (E3) is a retarded solution of the source-free wave equation, 
we see therefore that represents indeed a solution of the wave equation we had to solve : 

However this is not sufficient because we have also to solve the harmonic- 

(m) (m) 


coordinate condition (Eg). We shall refer to for the dehnition of an algorithm which 
permits to compute, simply from the algebraic and differential structure of the vacuum held 
equations, the necessary form of the second term in Eq. ([1.7| ), in such a way that the 
harmonic-coordinate condition be satished : = 0. In fact we shall not need, in the 

following, to be more precise about the latter term; simply we keep it into the form of a 
general retarded solution of the source-free wave equation, parametrized by some tensorial 
functions We assume that these functions are STF with respect to the multi-index 

L : i.e. so the multi-derivative Ol in Eq. i^X\) is a STF one (see Ref. 


for the notation). The latter construction represents the most general physical solution of 
the held equations outside the source . 

Let us now proceed with the formal re-summation of the post-Minkowskian series. That 
is, once the results (|2.5|) -( p.7| ) have been established for any order m, we sum them from 
m = 1 up to inhnity. In this way we obtain some formulas which are valid formally for 
the complete post-Minkowskian series, and, presumably, could hold true in a more rigorous 
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context of exact solutions. After summation we shall “forget” about the post-Minkowskian 
expansion, and consider that the exterior held represents merely the multipole decom¬ 
position of the actual held outside the compact support of the source (nevertheless, it is 
wise to keep in mind that the solution came from a formal post-Minkowskian summation). 
We denote the multipole decomposition by means of the calligraphic letter Af. Therefore, 
our definition is that the multipole expansion Af of the held outside the isolated source 
is nothing but the external solution constructed previously by means of the MPM method, 
and re-summed over the post-Minkowskian index m : 

MVn = >C, . ( 2 . 12 ) 

This dehnition is quite legitimate (and rather obvious) because we know that the MPM 
metric constitutes the most general solution for the exterior held. Thus, A1 (h^*^) is a solution 
of the vacuum held equations, now considered outside the physical domain of the source, 
r > a (while hg^t had been constructed for any r > 0). In that domain we have evidently 
the numerical equality 


(whenr>a). (2.13) 

After summation of Eqs. (|2.5|)-(|2.6|) over m, we get the near-zone structure 


VAeN, M(h^‘')= 5^hir“(lnr)PF£"„p(t)+ C>(r^) , (2.14) 

in which the functions A£0 p(t) = X]m=i ap(^)’ where a < N and p > 0. Notice 

that there is no lower bound for a because > —oo when m —>■ -|-cxo; similarly there 

is no upper bound for p. Secondly, coming to the constructive formula (p.7|) we obtain 

Mvn = \M{kn] + E 9 m (*-!:) . ( 2 , 15 ) 

;=o ' 
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where {t) = X]m=i following we shall regard the STF functions 

{t) as the “multipole moments” of the source, because they describe the physics of 
the source as seen from the exterior. We do not need to be more precise at this point. 
Let us simply comment that by imposing the harmonic-gauge condition (|1.2|) we hnd that 
there are only six components of these functions which are independent, and this yields 
the dehnition of six independent STF source multipole moments ■ ■ ■ (see Ref. 

m for the precise dehnition). Furthermore, the multipole-moment functions X^{t) have 
already been calculated in terms of the stress-energy tensor of a post-Newtonian source in 
Ref. m- However we prefer to leave these functions undetermined because we shall recover 
their expressions by means of a somewhat different method, and the agreement we shall hnd 
with the result of Ref. will constitute a crucial check of our computation. 


B. Near-zone expansion of the multipole decomposition 


In anticipation of the matching we consider next the inhnite near-zone re-expansion, 
when r —0, of the multipole expansion determined in Eq. (|2.15|) . We have already 

obtained the general structure of that expansion, given by Eq. (|2.14|) . Let us denote with the 
help of some overline the infinite near-zone expansion (without remainder), whose structure 
is therefore given by 


Wl(hA‘-)= ^hz.r“(lnr)^’F£;^(f), (2.16) 

where a G Z and p E N (and, of course, the multipolar index I eN). We must be careful at 
distinguishing the fully-hedged multipole decomposition which is dehned as soon 

as r > 0 and numerically agrees with the exact solution wherever r > a (in particular when 
r —>• -|-oo), from its formal near-zone re-expansion Al(h^^). Later we shall indicate the 
post-Newtonian expansion by means of the same overline notation. Indeed the near-zone 
expansion is really an expansion when r/A —0, which is equivalent to an expansion when 
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c —> + 00 , since the wavelength of waves is A = cP (with P = typical period of the internal 
motion). From the result (|2.15|) we can write 


= DRi [A1(A-'-)] + ^4 . (2.17) 

The overline in the second term means that one should expand the retardations t — r/c when 
r/c ^ 0. More explicitly we have 


d. 


j - r/c) 

I r 


^ (-)t . . ^ (i) 

i=o 


(2.18) 


where the superscript (j) indicates j successive time-derivations. The main problem is how 
to treat the hrst term in Eq. (p.l7|) . What we essentially want is to know how one can 
“commute” the operations of taking the near-zone expansion and of applying the retarded 
integral. In fact, the problem has already been solved in Ref. [^, which succeeded in 
writing the hrst term in Eq. (|2.17|) as the sum of an “instantaneous” operator, acting on the 
near-zone expansion of the source, and of a particular “anti-symmetric” wave (i.e. retarded 


minus advanced), solution of the source-free d’Alembertian equation. The result of Ref. 
Eq. (3.2), reads 


1-1 



4G^(-)'g 1 

^7^f(^-r/c)-7^^(t + r/c)l 



2r J 


(2.19) 

For completeness we present in Appendix ^ the proof of this result - a version of it which 
is somewhat improved with respect to that given in Ref. |^. The hrst term in Eq. (|2.19|) 
involves an operator acting on each of the individual terms of the formal near-zone 
expansion whose structure is given by Eq. ( |2.16|) , and which is essentially dehned by the 
solution of the wave equation that is obtained by iterated use of inverse Laplace operators, 
and regularized by means of our R-dependent hnite part procedure. Thus, 
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( 2 . 20 ) 


J-i 







where A ^ ^ = (A and the action of the inverse Laplacian on the generic term of Eq. 

( p.l6| ) follows from 


A-^ [hir^+“(lnr)P] = 




\dB) 

XB + a + 2- 1){B + a + 3 + l)_ 


( 2 . 21 ) 


[see also Eq. (|A16|) in Appendix The operator I ^ plays the central role in the present 
paper. It can be regarded as (the regnlarization of) the formal post-Newtonian expansion, 


when c —> -|-cxd, of the inverse d’Alembert operator, say = !/□ = 1/(A — ^d^). We 
can refer to as the operator of the instantaneons potentials, becanse it acts on the time 
variable t only throngh time-derivations, instead of involving a fnll integration like for the 
operator of the retarded potentials Notice that is closely related to the operator 

of the symmetric potentials, \ + '^Adv]! EH discnssion and the precise 

relation between these operators. As for the second term in Eq. (|2.19 ), it is made of an 
“anti-symmetric” wave, which represents in fact a solntion of the d’Alembertian eqnation 
that is regnlar in a neighbonrhood of the origin r = 0. Its near-zone expansion r/c —>■ 0 is 
composed only of terms containing some odd powers of 1/c : 


dj 


7 ^f(t-r/c)-7^f(t + r/c) 
2r 


H-C50 

-E 

i=l 


(2i+l) 


dL{r^^) 
(2i + l)! 


+ 00 


(2;+l)!! 




{2k+2l+l) 

„2k+2l+l 


( 2 . 22 ) 


k=0 


See also Eqs. (2.4)-(2.7) in Ref. [^ for alternative forms of the anti-symmetric wave. In 
the second of Eqs. (|2.22|) we have introdnced the nsefnl object 


A-’‘{xl) = 


( 2 / + 1 )!! 


(2A;)!!(2/ + 2fc + l)!! 


r^'^XL , 


(2.23) 
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which represents the iterated Laplacian operator A“^, regularized by means of the FP 


process, acting on xl which denotes the STF projection of the product xl = Xi^ 

[See also Eq. (|(Jllj|) for an alternative expression of the same object.] From Ref. [^, or 


Xil 

49 . 


, or 


from Eq. (|A11 ) in Appendix ^ we get the expression of the functions parametrizing the 
anti-symmetric waves, 


/ r-\-oo 

j dz^i{z) — z\x\/c) , (2.24) 

where |x| = |x|/ro [see Eq. (|2.9|) ]. These functions depend on the whole past-history of the 


source 


501. The 2 ;-integration involves the weighting function dehned by 




(2.25) 


This function is normalized so that dz ■ji^z) = 1, where the value of the integral is 
obtained by analytic continuation for I G C (see Appendix 0)- As shown in Ref. (see 
notably Section III.D there), the anti-symmetric waves in Eq. (|2.19| ) are associated with 
gravitational radiation reaction effects of non-linear origin. In particular they contain the 
contribution of wave tails in the radiation reaction force, which appears at the 1.5PN order 
relatively to the lowest-order radiation damping, i.e. 4PN order in the equations of motion 
m- To summarize this subsection, we have obtained the near-zone re-expansion of the 
multipole expansion as 


= J-i [>1(A/^^ 

+ 00 


4G (-)' a 1 Kr(i - r/c) - Kr(< + r/c) 

E—- 


1=0 


2r 


+ 00 


f - r/c) 


(2.26) 


1 r 

1=0 L 

The functions {t) are known from Eq. (|2.22|) , but the multipole-moments X^{t) have 
not yet been specihed at this stage (though they have already been calculated in Ref. [@). 
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Therefore, we have succeeded in computing the near-zone expansion as a functional 

of the sole unknown constituted by the X^^’s; only by matching can these functions be 
determined. 


III. INTERIOR FIELD 

A. Post-Newtonian expansion of the near-zone field 

Up to now, we have solved the Einstein held equations in the vacuum outside an isolated 
source (r > a), without any reference to the stress-energy tensor of the matter source. 
Our next task is to investigate the held equations inside and in the vicinity of the matter 
source, and more precisely in the so-called near-zone, or region for which r -C A, where 
A is the typical wavelength of the emitted waves. From now on we restrict attention to a 
post-Newtonian source, whose radius is a <C A. For post-Newtonian sources the near zone 
overlaps with the external region in what we shall refer to as the matching region, for which 
a < r -C A. In the matching region both the multipolar expansion of the exterior held and 
the post-Newtonian expansion of the inner held are legitimate. 

Let us denote by means of an overline the formal (inhnite) post-Newtonian expansion of 
the held inside the source’s near-zone : , which is of the form 

+ 00 .. 

V{x,t,c) = . (3.1) 

n=2 

By dehnition, the n-th post-Newtonian coefficient is the factor of the n-th power of 1/c; 

n 

however, we know from the structure of the near-zone expansion of the exterior held [see 
Eq. (|2.16|) ] that the post-Newtonian expansion will involve also, besides the usual powers of 
1/c, some logarithms of c (in fact when stating this we are anticipating on the result of the 
matching). So the coefficients still depend on c via the logarithm of Inc, and from Eq. 

n 

( ^.l(j| ) we infer that they are in fact some power series in In c. The hrst appearance of In c is 
at the 4PN order (i.e. corresponding to a term ~ Inc/c® in the equations of motion) and is 
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associated to the physical effect of wave tails |^. In Eq. (|3.1|) we have indicated that the 
expansion starts at the level 1/c^, but we could be more precise because the Oi component 
of starts only at the level 1/c^, while the ij component is at least of order 1/c^. This 


does not matter for our purpose; simply in our iteration we include these post-Newtonian 
coefficients as zero : = 0, and Ti^ = TT 

2 2 3 

( lOI) we have the same type of expansion, 


coefficients as zero : h =0, and h = h =0. For the total stress-energy pseudo-tensor 


+ 00 




n=—2 

r.2 


(3.2) 


The expansion starts with a term of order corresponding to the rest mass-energy of the 
source has the dimension of an energy density). Here we shall always understand the 
infinite sums such as ( p.l|) - (|3.2|) in the sense of formal power series, i.e. merely as an ordered 
collection of coefficients : e.g. () . We do not attempt to control the mathematical 

nature of these series. 

In this article we make two important assumptions. First, we assume that the post- 
Newtonian coefficients (and similarly are smooth functions of space-time. 


e . (3.3) 

n 

Evidently this comes from our consideration of regular (smooth) extended matter distri¬ 
butions, described by T^^(x, t) G C“(M^), a priori excluding black holes or point-particle 
singularities. Second, we assume that the structure of the expansion at spatial infinity, i.e. 
r —>• -|-cx) with t = const, is of the type 


VAfeN, + 


(3.4) 


(and similarly for each . On purpose we have written an expansion which is very similar 

n 

to the one in Eq. (|2.14|) , because as we shall see the functions Fi^n,a,pit) will be equal to 
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the post-Newtonian coefficients of the functions FL^a,p{t) appearing in Eq. (p.l4|) . However 
it is important to realize that in contrast to Eq. (|2.14|) which is a near-zone expansion [cf. 
the remainder (9(r'^)], the expansion written in Eq. (|3.4| ) is a far-zone one, with remainder 
0{l/r^). It would have been clearer to write the latter expansion with some (lnr)^/r^ 
with b = —a, but since we are going to show, from the method of matched asymptotic 
expansions, that the inhnite far-zone expansion (ignoring the remainder) is actually the 
same as the inhnite near-zone expansion, it is better to write it in this form, with the range 
of the powers of r in Eq. (|3.4|) being —a < N instead of a < in Eq. (|2.14|) . In doing so we 
are again anticipating on the result of the matching. Finally, we assume that, at any given 


post-Newtonian order n, the maximal divergency of the far-zone expansion (|3.4|) is hnite, 
i.e. there exists some aniax(^) £ N such that a < amax(?^)- 


Next we perform the iteration of the post-Newtonian held (|3.1| ) up to any order. Our 
strategy consists of hnding the general post-Newtonian solution of the relaxed Einstein held 


equation ( p..l|) . This solution will depend on some arbitrary “homogeneous” solutions, in the 
form of harmonic solutions solving the source-free d’Alembertian equation (in a perturbative 
post-Newtonian sense). In a second stage we shall obtain these harmonic solutions by 
imposing the matching to the external multipolar held obtained in Section Finally we 
shall check that our post-Newtonian solution is divergenceless, i.e. it satishes the harmonic- 


coordinate condition (|1.2|) , in consequence of the conservation of the stress-energy pseudo 
tensor . Notice that we do not try to incorporate into the post-Newtonian series the 
boundary conditions at inhnity [viz the no-incoming radiation condition). Indeed this is 
impossible at the level of the post-Newtonian expansion considered alone, because its validity 
is limited to the near-zone. Even if we dehne an “improved” post-Newtonian series by 
considering some retarded integrals that are formally expanded when c —^ -|-cxd as series of 
Poisson-like integrals |^, we ultimately end up with an inconsistency, because the Poisson- 
like integrals are some local-in-time functionals, depending on the source only at the current 
time t, and we know that the post-Newtonian held starts to depend on the whole past 


history of the source from the 4PN order ||31|-p3|. Therefore, we do not follow this route in 
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the present paper, and, instead, we incorporate into the post-Newtonian series the boundary 
conditions concerning the wave held at inhnity by means of the matching equation. 


We insert the post-Newtonian ansatz (|3.1j)-(p^ into the “relaxed” Einstein held equation 


(p..l|), and equate together the powers of 1/c. The result is an inhnite set of Poisson-type 
equations : 


Vn > 2 , Ah'*'' = . 

n n—4 n—2 


(3.5) 


Evidently, the second term comes from the split of the d’Alembertian operator into a Lapla- 
cian and a second time derivative : □ = A — the time derivative c?o = is smaller 
than the spatial gradient di by a factor 1/c - this is the basic tenet of the approximation. 
When n = 2 and n = 3 the second term in Eq. (|3.5|) is zero, which we take into account by 
assuming that = 0. We proceed by induction, i.e. we hx some post-Newtonian 

order n, assume that we succeeded in constructing the sequence of previous coefficients 
, • • • , and from this we infer the next-order coefficient . 

2 n—1 n 


The most general solution consists of the sum of a particular solution and of the most 
general admissible solution of the homogeneous equation, which is simply the source-free 
Laplace equation. Let us hrst find a particular solution. We recalled in the introduction 
that the usual Poisson integral cannot be used to dehne a solution, because the bound at 
inhnity becomes rapidly divergent when going to higher and higher post-Newtonian orders. 
Fortunately, thanks to our two assumptions and 0). we shall be able to dehne 

a generalized notion of Poisson integral, in a way similar to our previous dehnition of a 
retarded integral operator in Eq. (|2.10|) . That generalized Poisson integral will constitute 
an appropriate solution of the post-Newtonian equation. For any source term like which 

n 


is at once smooth, Eq. (|3.3|), and admits a far-zone expansion of the type ( |3.4|) [note that 
Eqs. ( |3.3|) - (|3^ ) hold for as well as for we multiply it by the same regularization 
factor as in Eq. ([^.lUD, and then apply the standard Poisson integral. The result. 
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(3.6) 


A-‘ 


n 


(x,t) 


1 

Att 



<Py 

x-y 


n 




where lyl*® = ly/rol^, defines a certain fnnction of i? G C. The well-dehniteness of that 
integral heavily relies on the behavionr at the bonnd at inhnity. There is no problem with 
the vicinity of the origin becanse of the smoothness of the integrand. From the asymptotic 
expansion (|3.4|) , with a < Omax [recall that Omax = amax(^)], we hnd that the integral 
converges at inhnity when < —Omax — 2. Next we can prove that the latter fnnction of 

B generates a (nniqne) analytic continnation down to a neighbonrhood of the origin B = 0, 
except at i? = 0 itself, at which valne it admits a Lanrent expansion with mnltiple poles np 
to some hnite order. More details are given in Appendix p. Then, we consider the Lanrent 
expansion of that fnnction when B ^ 0 and pick np the hnite part, or coefficient of the 
zero-th power of B, of that expansion. This dehnes onr generalized Poisson integral : 


^ FP 

L n B=0 




(3.7) 


The hnite-part symbol FP has exactly the same meaning as in Eq. (p^.lU| ). However, notice 
that in contrast to Eq. (|2.10|) where the regnlarization factor dealt with the singnlarity 
when r —> 0, and hence snpposes initially that ^{B) is a large positive nnmber, in Eq. (p.7|) 
the regnlarization concerns the behavionr of the integral when r —> +cxd, and so one mnst 
start with the sitnation where ^{B) is a large negative nnmber. The main properties of onr 
generalized Poisson operator is that it solves the Poisson eqnation. 


(3.8) 

and that the solntion owns the same properties as the sonrce i.e. the smoothness, 

n n 

Eq. ( |3.3|) , and the particnlar far-zone expansion given by Eq. (|3i^). These facts are proved 
in Appendix Therefore, we have fonnd a particular solntion of the Poisson eqnation, and. 


A 


A-i 
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furthermore, this solution can be iterated at will, because the operator keeps the same 


properties from the source to the corresponding solution. Quite naturally we denote the 


iterated Poisson operator = ( A“^ 


fc+i 


; it is not difficult to show that 


A-'=-i[rn(x,t) =-—FP / d^y 
n 47r B=0 Jm 3 


|X- y 


2k-l 


(24)! 




(3,9) 


From that integral we obtain the operator of the “instantaneous” potentials exactly in the 
same way as in Eq. (|2.20|), but now acting on post-Newtonian coefficients such as i.e. 

n 

satisfying both Eqs. ( p.3|) and (|3.4|) : 


+ 00 


p\ \2k ^ __. 


fc =0 


(3.10) 


It is clear that we have a particular solution of d’Alembert’s equation : 


□ 




n 


— 


(3.11) 


We can check that the dehnition we have proposed in Eqs. (p.20|) -( ^ATD is a particular case 
of the more general dehnition (p.9|) - (|3.10|) . Indeed, if we apply the formulas (|3.9|) - (|3.10|) to 
one of the terms composing the “far-zone” expansion of the post-Newtonian coefficient, i.e. 
hLr“(lnr)^F’(t), we get the same result as the one resulting from Eqs. (|2.20|) - (|2.21|) . 

By means of the Poisson operator A“^ so constructed we hrst hnd a particular solution 
of Eq. 


(= leTrG . (3.12) 

/ part n—4 n—2 

To this solution we add the most general solution of the homogeneous Laplace equation. It 
can be written, using the STF language, as the sum of two multipolar series, one of them 
being of the type xl, that is regular at the origin r = 0, the other one being like dii^/r), 
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i.e. regular “at infinity” r —>■ +cxo (see Ref. for the notation). Imposing the smoothness 
condition ( |3.3|) for the post-Newtonian held we discard the second type ~ diiX/r), and 
retain as the only admissible homogeneous solution the hrst type ~ xl- Therefore, we hnd 
that there must exist some STF-tensorial functions of time, say such that 


+ 00 

=^B^j^{t)xL . (3.13) 

^hom ^ n 

The functions B^{t) will be associated with the reaction of the held onto the source, and 

n 

depend on which boundary conditions are to be imposed on the gravitational held at inhnity. 
The most general solution for the n-th post-Newtonian coefficient thus reads 


r'' = (r-) 


part 


" (r) 


horn 


(3.14) 


It is now trivial to iterate the process. We substitute for in the right-hand-side of Eq. 


n—2 


( 3 .I 2 I ) the same expression but with n replaced by u — 2, and similarly descend untill we 


stop at either one of the coefficients = 0 or = 0. At this point W is expressed in 

0 1 n 




-flU . 


terms of the “previous” and B^'s, with m < u — 2, i.e. 

m 


[f]-i 


_ _ +00 [ 2 I 1 (2k) _ 

/r = 161G ^ j + Y.Y. ■ (3.16) 

n ^^ n-A-2k ^^ ^^ 

fc=0 1=0 k=0 

Here [n/2] denotes the integer part of n/2; means the 2A;-th partial time derivative 
{d/dtY^ and the superscript {2k) the 2fc-th total time derivative; the operator is the 

one dehned by Eq. (|3.9|) ; the object A~’^{xl) has already been introduced in Eq. 


Once we have the result (EH)’ we “re-sum” it from n = 2 up to inhnity. After commuting 
the summations over n and k we arrive at 


^P.-rrC- _ 1 {2k) _ _ 

r- = + Y.I: . 




1=0 k=0 


(3.16) 
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where we have recognized the operator of the “instantaneous” potentials as dehned by Eq. 
( p.lOj ), and where the functions read 


+ 00 


A more compact alternative form is 


( 3 . 17 ) 


n=2 


fc"" = PI + E ^ 21-1 |Br(i) iil 

1=0 


( 3 . 18 ) 


Actually the latter forms are not the best for our purpose. Since the hrst term in Eqs. 
( p.l(j| ) or ( |3.18|) is a particular solution of the d’Alembert equation [see Eq. ( p.ll| )], the 
second term is necessarily equal to (the near-zone re-expansion of) a homogeneous solution 
of the source-free wave equation, and most importantly a regular solution at it. So it should 
be in the form of some anti-symmetric multipolar waves : retarded minus advanced. Indeed, 


this readily follows from the second equality in Eq. (|2.22|) . We introduce a new dehnition 
^r(^) by posing 


BTit) = 


( 21 + 1 ) 

Af (t) 


(3.19) 


c2^+i(2/ + l)!! ’ 

where the /-dependent factor is chosen to match with Eq. (|2.22|) . [Because of our assumption 
of stationarity in the past, t < —T, the relation ( p.l9|) determines (t) up to a constant. 
However it is clear that this constant will cancel out in the anti-symmetric wave in Eq. 
( p.20| ).] In terms of this dehnition, we hud the hnal result of this section. 




H-OO 


+5^4 


[ <(i-r/c)-3lf(i + r/c) 

2r 


(3.20) 


1=0 I 

where we recall that the overline means the post-Newtonian or equivalently near-zone expan¬ 
sion [see Eq. ( |2.22|) ]. For the time being we shall refer to the {t)A as the radiation-reaction 
functions. 
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B. Multipole expansion of the post-Newtonian solution 


In the previous section we obtained the general solution for the post-Newtonian expansion 
in the form (p.20|) , and parametrized by some (for the moment) unknown radiation-reaction 
functions A'^{t). To arrive at this we made an assumption concerning the particular struc¬ 
ture for the far-zone expansion, at spatial inhnity, of the post-Newtonian coefficients : Eq. 
O' Here we shall denote the corresponding inhnite expansion (without remainder term) 
by means of the same calligraphic letter M. as used to denote the multipole expansion, be¬ 
cause the far-zone expansion of the post-Newtonian coefficients is equivalent to a multipolar 
decomposition. From Eq. ( 3.4|) we have 


M{h^'')= 5^hir“(lnr)PF£‘; (f). (3.21) 

So, summing up the post-Newtonian series. 




(3.22) 


where the functions involved are Fi^a,pit) = '^n =2 ^PL,n,a,p{f)■ As we can see, the far-zone 
expansion that we have just postulated is exactly the same, with the same functions FL^a,pit), 
as the near-zone expansion we had previously written in Eq. (p.l6|) . This equality is already 


the matching equation between the near-zone expansion of the multipolar held, Af (h^^), and 


the multipolar-far-zone expansion of the post-Newtonian held, A4{h^’^), whose consequences 


will be investigated in Section [TV. 

The fundamental result which is needed for computing the far-zone expansion of the post- 
Newtonian series concerns the expansion of the generalized integral operator acting on 
the post-Newtonian source More precisely, we are interested in knowing under which 
conditions one can commute with the operation Ai of taking the far-zone expansion. 
Clearly, the two operations can be commuted at the price of adding some homogeneous 
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solution of the d’Alembert equation. We prove in Appendix that the latter homogeneous 
solution is made of multipolar waves of the symmetric type, i.e. retarded plus advanced. We 
obtain 


__ -| +00 / W 

M (l-i [t""]) = I-i|A<(r'‘-)] - T- E 

^ i=0 

Here the overline notation has the same meaning as in Section |1| : this is the Taylor expansion 
when r —^ 0, but that expansion should be considered, a posteriori, as an expansion when 
r —>■ + 00 . That is, our notation means 


.Fr(t-r/c)+JT(t + r/c) 


2r 


(3.23) 


dj 


.Fr(t-r/c)+.Fr(t + r/c) 

2r 


+ 00 

E 


(2i) 




{2i)\ 6 

dL{r~^) 


2 = 0 
-\-oo 


{2i) 


r,2i 


(3.24) 


i=0 


[see also Eq. (|C14|) ] where the right-hand-sides are to be considered as some expansions at 
spatial inhnity, of the general type given by Eq. ( 3.22|) . The functions {t) parametrizing 
these symmetric waves are STF and explicitly given by 


n^it) = 



j |5|^A-t [xl] t) , 


(3.25) 


where A”-? [xi\ is given by Eq. ( |2.23|) . See the proof in Appendix y. An alternative form 
reads as 


= FP J d^x|x|'®h 2 , J dz 6i{z) T^''{x,t — z\x.\/c) , (3.26) 

where the integration over the z-dependent cone t — - 2 |y|/c involves a weighting function 
6i{z) that is closely related to the function 'yi{z) introduced in Eq. (|2.25|) : 


Hz) 


(2/ + 1)!! 

2i+H\ 


(l-zH 



(3.27) 


28 






















and whose integral is normalized to one : dz6i{z) = 10. The fnnction 6i{z) approaches 

the Dirac delta-fnnction in the limit of large I : lim i^+oodi{z) = S{z). In Eq. (|3.26|) we 
have indicated by means of an overline the fact that this expression is valid only in a sense 
of post-Newtonian expansion. Note that becanse the latter post-Newtonian expansion is 
“even”, containing only even powers of 1 /c, one can replace the argnment t — ^|y|/c inside 
eqnivalently hj t + z\y\/c. 

Finally, thanks to Eqs. (|3.23|) - (|3.27|) , we are in a position to write the inhnite multipolar- 
far-zone expansion of the post-Newtonian solntion as 

IOttG; 






1=0 


2r 


+ 00 

1=0 


I klf (t-r/c) -Af (t + r/c) 

I 2 r 


(3.28) 


We recall that the radiation-reaction fnnctions (t) are still nndetermined at this stage. 
The symmetric and anti-symmetric waves are given by Eqs. (p.24|) and (p.22|) respectively, 
considered here as inhnite far-zone expansions. 


IV. MATCHING 


In Section EI3, we fonnd the most general expression for the mnltipolar expan¬ 
sion satisfying the no-incoming radiation condition, in terms of some nnknown 

“multipole-moment” STF-functions {t) [see Eq. ( p.l5|) ]. On the other hand, in Sec¬ 
tion |111 A|, we obtained the most general solution for the post-Newtonian expansion 


^The normalization for the function 6i{z) is consistent with that of the function 7 ^( 2 ) : 

dz ji(z) = 1 , owing to the fact that the integral f^^dz (1 — z'^Y is zero by complex an¬ 
alytic continuation in / G C. 
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as parametrized by a set of unknown “radiation-reaction” STF-functions (t) [see Eq. 
(p.20|)]. We are now imposing the matching condition 


=M(h‘ 




(4.1) 


In fact we have already postulated this equation when writing that the two formal expansions 
( p.l6| ) and (p.22|) are the same. Recall that the matching equation ([4.1|) results from the 
numerical equality verified in the exterior near-zone : a < r X. It is 

physically justified only for post-Newtonian sources, for which the exterior near-zone exists. 
The matching equation is actually a functional identity, i.e. true V(x, t) G M* x M; it 
identihes, term-by-term, two asymptotic singular expansions, each of them being formally 
taken outside its own domain of validity. In the present context, the matching equation 
insists that the infinite near-zone expansion, r —0, of the exterior multipolar field is 
identical to the inhnite far-zone expansion, r —> -|-cxo, of the inner post-Newtonian held. Let 
us show now that Eq. (|4.1|) permits determining all the unknowns of the problem : i.e., at 
once, the multipole moments Xff' and the radiation-reaction functions . In particular 
we hnd that the multipole moments Xff are in agreement with the earlier result derived in 


Ref. \i0 


For the sake of clarity we re-state here the two results we reached for the two sides of 
Eq. (^T|). The left-hand-side was obtained in Eq. (p.26|) : 


~ l\ 

1=0 




2 r 


+ 00 


f AT(t - r/c) 


(4.2) 


1=0 1 

in which the functions TZf', which come from the non-linearities of the held equations in 
vacuum, are known from Eq. (|2.24|) . The right-hand-side of the matching equation was 
found in Eq. (p.28|) : 
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Mih^") = 


IGttG 




^ ^~[r 

1=0 


dr 


J^^,\t-r/c)+T^\t + T/c) 
2r 


+ 00 


/=0 


[ Af (t-r/c) - v4f (t + r/c) 


E®. , 2 , 


( 4 . 3 ) 


Here, the functions which depend on the matter and gravitational content of the post- 
Newtonian source, take the dehnite expression given by Eqs. (|3.25|) - (|3.26|) . 

Comparing the equations ( |4.2|) and ( |4.3| ), we readily discover that they share an obvious 
common term, that is the hrst one. Indeed, we manifestly have 


X 


-1 


2W(A 


fiiy] 




( 4 . 4 ) 


The hrst equality comes from the matching equation, as applied to the gravitational source 
term and the second equality comes from the fact that the matter tensor has a 


compact support, so that = 0. Hence the two hrst terms in Eqs. (^]2|) and (|4.3|) 

match together. This is a somewhat remarkable fact, because most of the complexity of the 


Einstein held equations is actually contained into these terms, either ] for the 

external held or for the inner one. But for doing the matching we don’t 

need all this complexity; these two terms match and therefore are to be identihed. Notice 
also that this is a non-trivial result, since the two sides of Eq. (|?T) strongly depend on the 
yet unknown functions and , which enter the latter two terms in a very intricate 
way, coupled together as they are via the non-linearities of the held equations. Nevertheless, 
the matching equation tells us that these terms must be rigorously identical. 

As soon as we have noticed that the hrst terms in Eqs. ( |4.2[) and ( [4.3|) are equal, we can 
compare the other ones, and because the retarded and advanced waves have some diherent 
structures, they must be matched independently, so we get two relations to be satished. We 
hnd that these are solved if and only if the multipole-moments in the exterior held and the 
radiation-reaction functions in the inner held are given by 
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(4.5) 


xr(t) = 

AT it) = 


4G'(-) 


c4 l\ 
4G(- 


-tFTit) , 


r 


l\ 


Ffit) +Kr(i) 


(4.6) 


Therefore, both the multipole-moments and the radiation-reaction terms are determined as 
some explicit functionals of the pseudo-tensor and nothing else. (Actually we could add 
any constant to the definition of AT (t), but this is physically irrelevant because the constant 
disappears from the anti-symmetric waves; see also Ref. p0| .) 

Finally, by way of summary of the results, we take back the latter expressions and fill in 
the external and inner fields, which are then entirely determined as coming from a unique 
solution of the Einstein field equations in harmonic coordinates, valid everywhere inside and 
outside the source. The exterior field is 


4G 


+ 00 


iA4(A'“')i - ^ E IT®' 

1=0 


n-O - r/c) 


(4.7) 


where the multipole moments are given in terms of the post-Newtonian expansion of the 
stress-energy pseudo-tensor by 


.Fr(t) = FP J d^ylyTvL j^^dz6^{z)T>^Ty,t-z\y\/c) 

+ 0O .. „ 

= E I lyi-] ‘) ■ 

A —n 


(4.8) 


This result is in perfect agreement with the multipole decomposition of the exterior field 
obtained in Ref. [see Eqs. (3.13)-(3.14) there]. On the other hand, the inner post- 
Newtonian field is given by 




4G g (-)' f ATit - r/c) - AZit + r/c) 


1=0 


2r 


(4.9) 


where the radiation-reaction function is composed of two terms : 


AT{t) =.FD*)+Kr((). 


( 4 . 10 ) 
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The first term is nothing but the exterior multipole moment given by Eq. ([4.8|) , and one 
can check that it contains the standard radiation-reaction effect at the 2.5PN order. The 
7?.^^-term is dehned by Eq. (|2.24|), or, rather, the post-Newtonian expansion of it, i.e. 


TiTit) = FP J Sy |y 


B 


b + CXD 




( 4 . 11 ) 


This term is quite interesting : it depends on the non-linearities of the exterior held, de¬ 
scribed by the gravitational source term (or, more precisely, the non-stationary part 

of it IQ), which are to be computed by means of the multipolar-post-Minkowskian algorithm 
of Refs. p6|j3^ (see in particular Section III.D in Ref. [Q for some detailed computations 
of this term). Physically the function contains the effect of wave tails in the radiation 
reaction force which arises at the 4PN order [^Tj-Q. It is not difficult [using notably the 
formula (^.21|) below] to derive the more explicit expression for the contribution of to 
the anti-symmetric wave in Eq. 


TZTit + rM l 

j^ira^FP /<lV|yr|^9fP-‘>f(r-)(y.i-|y|/c), (4,12) 

When they are computed by post-Minkowskian approximations, the remaining integrals 
will typically yield, after integration over the angles, some “hereditary-like” contributions, 
depending on the whole integrated past of the matter source (see Ref. [Q). 

It is tempting to speculate that the second term in Eq. ( |4.9|) , made of the anti-symmetric 
multipolar waves parametrized by the functions can be regarded as the contribution, 

in a sense to be made more precise, of the radiation reaction forces at work inside the 
post-Newtonian source. [Indeed we have checked that these functions contain the known 
radiation-reaction terms at the dominant 2.5PN order as well as the dominant contribution 
of tails at the 4PN order.] We shall leave for future work the systematic study of this term, 
as well as the possibility to answer the latter speculation. 


dr 


nnt-r/c)- 


E 


Oii\(l - 441 


2 r 

+ 0 
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V. HARMONIC-COORDINATE CONDITION 


The latter solution for the post-Newtonian held, Eqs. (|4.9|) - (|4.12 ), has been obtained 
without imposing, in an explicit way, the condition of harmonic coordinates (|1.2|) . Indeed, 
we have assumed this condition to be true, and we simply matched together the post- 
Newtonian and multipolar-post-Minkowskian expansions, satisfying the relaxed Einstein 
held equations (0) in their respective domains. We found that the matching determines 
uniquely the expressions of the multipole moments {t) and radiation-reaction functions 
(t) as some functionals of the stress-energy pseudo-tensor However, we never used 
the harmonic-coordinate condition during the matching; it was not necessary for the formal 
determination of the unknown parameters {X^, A'^). Therefore, it is quite important to 
check that our post-Newtonian solution is divergenceless as a consequence of the conservation 
of the pseudo-tensor [see Eq. (EDI. SO that we really grasp a solution of the full Einstein 
held equations. 

We check the divergenceless of directly on Eq. (p.l 6 |) . We apply the operator on 
each side of the equality : 


Bjr = [t'“’1 + a. 


r +O0 


1 m _ 


^2k 


( 6 . 1 ) 


We must transform the two terms on the right-hand-side in order to make explicit the fact 
that these two terms are exactly the opposite. The hrst term, that is to say, the divergence 
of the operator, is not obvious since, even if time derivatives commute with spatial 
derivatives do not. 


= 2:-' [dof^^] + [r^] . (5.2) 

is a sum of and spatial derivatives do not commute with because of 

the |y|^ factor [see Eq. (|3.9|) for the exact expression]. To see how to tackle this problem, 
let us start with the spatial divergence of We assume that r(x, t) is a function of the 
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“post-Newtonian” type, i.e. satisfies the requirements (|3.3|) - (|3.4|) . The di derivative, in Eq. 
( p.3|) , applies hrst to the fc**^-Poisson’s kernel but after having noticed that the x*-derivative 
of this kernel was equal to minus the ?/*-derivative of it, we can make an integration by part 
and distribute the ?/*-derivative on |y|^ and on r so that 


= -Tf_p 


2fc-l 


7-(y,t) = 


47riJ=o y {2k)\ 


2k-l 




2k-l 


-r(y,t). (5.3) 


The hrst term in the last line of the previous equality is equal to A~^~^{dif). Now, let 
us concentrate on the last term in the same line. We can, of course, write di{\y\^) = 
Bni\y\^~^/ tq. Moreover, since T{y,t) is regular at the origin (|y| = 0), the integral is 
always convergent on any neighbourhood of the origin. Translating these two remarks in 
the last integral of Eq. ( |5.3|) , and since we take the hnite part when B = 0, this last integral 
is zero, because of the explicit factor B, when ranging from |y| = 0 up to some arbitrary 
hnite value |y| = TZ. So, after replacing |x — yp^“^ by {2k)\ A“^(|x — y|“^) we are left, in 
Eq. (|5.4|) , with one integral ranging over |y| >71, 

I2k-1 


FP 

B=0 ^ 

= FP 

B=0 


rfV.dyl 


S' 


|x-y| 


i2k)\ 




I d-^yOmnM A-^(|x 

|y|>7^ 




V Vi^A^-(a;i)FP 

l\ '' ^B=0 


M{T){y,t) 

dV5*(|y|^)^"0L(|y|-'))-M(r)(y,t). 


(5.4) 


n=0 l>0 "■ ■'\y\>'^ 

In the last line of the previous equation, we expanded the fc*'^-Poisson’s kernel for |y| 3> |x| 
using Eq. ( P12|) . This is possible thanks to the fact that 71 is arbitrary and may be chosen 
such that TZ >> |x|. We also note that r turned into M.{t) because r = A4 (t) in the far 
zone. In this way. 


diA-'^-’^ [7‘(x,f)] = A-i-^ [diT{x,t)] 


tAAK'"(ii)FP / dVa.(|y|«)A->(Si(|yr‘))M(r)(y,(). (5-5) 

^ n=0 l>0 ■ dlylyn 
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and we notice that the commntation of the spatial derivative and the generalized /c*^-Poisson 
integral depends only on the behavionr of r(x, t) at spatial inhnity. This fact was foreseeable 
since for a fnnction r(x, t) with compact snpport the commntation wonld be trivial. Thanks 
to the general resnlt given by Eq. (|5.5|) , in which we replace r by we can determine the 
spatial divergence of We can then get d^I~^ that is the snm of 

and a non-trivial term. Since = 0, the resnlt for the hrst term of Eq. (|5.1|) rednces to 

the non-trivial term, that is to say : 


It*-] 

r.4 /I ^ B=0 


(2k+2n) 


C’ — l\ 

k,l,n>0 


|y|>^ 


d^ydmnA-HdLi\y\-^))Mir^){y,t) . (5.6) 


Now, we want to prove that the second term in the right-hand-side of Eq. ( |5.1|) is exactly 
the opposite. In Eq. (p.7|) , we expand this last term in its 3+1 form so that we can treat 
separately terms with time derivative and terms with spatial derivative. 


1 (2«) „ 

n ,/>0 


( 2 n ) 


[XL 


+ E 


{XL 


( 6 . 7 ) 


n ,/>0 


The hrst term of Eq. (E3). thanks to a STF formnla, can be written withont the nse of 
spatial derivative. The index i coming from this derivative is distribnted on the mnlti-index 
L in the way 


1 ( 27 *) „ 


n ,/>0 




n ,/>0 


2/+ 3 


( 2 n ) ( 2 n ) _ 


^{xL-l) 


(5.8) 


In the second term of Eq. (|5.7|) , we express the fnnction in terms of and [cf. 


Eqs. ( |3.19|) and (|4.6| )] becanse the time derivative, do, will act on the integrand of these two 
time-varying moments : 


_ AG(-Y — r 'i 

A-H£L)doB^L^{t) = 

First, we investigate the case of doB^, nsing the formnla (p.25| ), where the time derivative 
acts on , 
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( 5 . 10 ) 


I2li) 


A-"(fi)5„.F““(() = A-"(fi)5^^FP / <iVly|®A-'-(fc)S„T»''(y,() . 


k>0 


We can replace by —diT^’^ thanks to the conservation eqnation of the psendo-tensor. 

After integrating by part we get 

(2k) 




fe >0 


(2k) 


+A-”(ii) ^ ^FP / <iV|y|®Si(A-‘(yi))y‘'(y.i). 


n >0 


( 6 . 11 ) 


The same STF formnla as nsed in Eq. (^) enables one to transform the second term of 

Eq. (|5.11|) so that, at the end, we get the dehnitive resnlt : 

_ _ _ 1 r _ _ (2fc) 

= A-"(hL)^—FP^ / d?ydi{\yf)^-^{yL)T^'^{y,t) 
k>0 ^ 


.- A“^frri 


(5.12) 


c2(2/ + 3)' 

We can, now, investigate the case of the hrst term in Eq. (|5.9|) , which is a little bit more 
complicated since it involves a retarded integral. 


A-"(hL)5o7^i (t) = A-^(hL)|T y d y |y| y dz'^i{z)M{doT^''){y,t-z\y\/c), (5.13) 

where the fnnction 7 z(^) is given by Eq. (|2.25|) (for simplicity’s sake we do not write the 
overline indicating the post-Newtonian expansion). We do the replacement of into 

—diT^^. Before integrating by part, we better have to notice that the partial derivative di 
acts on which is then evalnated at the event (y,f — 2 ;|y|/c); we mnst be carefnl abont 
the space dependence of the time variable t — z\y\/c. The last eqnation then becomes 

-A-^(hL)FP [ d^y\y\^yLdi( [ dz-fi{z)M{T^''){y,t - z\y\/c) 


5=0 


( 1 ) 


A-'^(hi)FP / dV|y| yirii dz-'yi{z)M{T^’'){y,t-z\y\/c). 


5=0 


( 6 . 14 ) 


In this way, the hrst term can be integrated by part straigthforwardly, in terms of d^y 
integration, showing np a di{\y\^) term and a di{yL) term. The second term will also be 
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integrated by part, in terms of dz integration, using the fact ^[ 7 /+i( 2 ;)] = —{21 + 3)z'yi{z); 
so we have 


A-’"(£i)|T J d^ydi{\y\^)yL J dz-fi{z)M{T^’'){y,t - z\y\/c) 
+A-"(fL)FP [ d^y\y\^di{yL) [ dz-fi{z)M{T^’'){y,t - z\y\/c) 

B=0 I /i 


A-^(£l) 


FP / d^y\y\^yLU, 


c^{2l + 3 ) 5=0 


( 2 ) 

*71+1 WX(T“')(y, t - z\y\lc) 


( 5 . 15 ) 


The sum of these three terms can be transformed so that the function TV'j^ shows up. Since 
for any STF tensor fLdi{yL) = ifn-iyL-i and fiyLyi = fiyiL + i^fiL-iyL-i\y\^, and 
keeping in mind that all the multi-indices L will have to be summed, we can write 


A-’^(xl)FP / dydi{\y\ )yL / dzyi{z)M{T^’'){y,t - z\y\/c) 

5=0 / /i 


(5.16) 


+l A-'^{xiL-i)FF^ J dy\y\^yL-i 

, ^^{xl) 


dzyi{z)M{T^'"){y,t- z\y\/c) 


( 2 ) 


^ 2 ( 2 / + 3 )'|Jo j / dzyi+i{z)M{T^^){y,t-z\y\/c) 




( 2 ) 


FP / dV|yr|yr^L-i / dzyi+,{z)M{rn{y,t-z\y\/c) 


c2(2/ + l)(2/ +3)5=0 

An interesting relation between 7 rfunctions : -£j[''yi+i{z)] = (2/-|-1)(2Z-|-3) [ 7 i_i(^) — 7 /( 2 ;)], 
after integrating by part the last integral, in terms of dz integration, allows us to get the 
more explicit form : 


A-^{xL)donl^{t) = A-^{xl)FF / d^ydi{\yf)yL / dzyi{z)M{T^‘'){y,t-z\y\/c) 

B=0 _/ _/i 


-- A~'^(xr) 

+/A-(£,^_l)7^^_l(^) + X it) • 


c2(2/ + 3) 


Summing up Eqs. (^. 121 ) and (|5.17|) we obtain 


,__ (2n) __ ^ 97 —V __ 

«>ci l,n>0 

(2n-|-2Z-|-l) 

+ FP f <iVa.(|y|®)»L /'”*7 i(2)'^*^'"* 


Eg / <'VSi(|y|®)A-”(»i) 


k>0 


( 5 . 17 ) 


(2n-|-2fc+2Z+l) 

^2n+2fc+2Z+l 


iy,t- z\y\/c) 


y — 


^2n+2/+l 

-1 {2n) {2^) 


n,/>0 


2 /+ 3 


(5.18) 


38 















The last line cancels ont the terms coming from Eq. (^.81) . 

We can therefore write down the result for the divergence of which, at this stage, 
depends only on terms with integrals of di{\y\^) and having the spatial structure given by 
After summing Eqs. (^) and (|5.18|) we get 


Of^h = 

4G y. 
~~ 2-^ 


(5.19) 




FP 

C’ —' l\ c^k+2n 

n,/,/c>0 


d'ydi{\y\^)A-f^{dL{\y\ ^)) Al(r*^)(y,t) 


(2n+2fc) 


'\y\>n 


\i 


A-^{xl) 


4G sy- 

fA Z_^ + 1)! C^"-+2fc+2*+l B=0 y|v|>77 

n,l,k>0 ^ 


FP 


d^yd,{\y\'')^-KyL)M{rn{y.t) 


{2n+2k+2l+l) 


4G 




\i 


A-^{xl) 


{2n+2l+l) 


C’ —' (21 + 1)! B=0 

n,l>0 ^ ' 


FP / d^ydmnvL / d^7z(^)A^(0(y,^-^|y|/c) • 


In the second term we have used the fact that the integral depends only on the values for 
which |y| > 7^ to write on that domain. The last term of Eq. ( b.l9|) depends 

on a retarded integral of the multipolar post-Minkowskian expansion By integrating 

by part the integral over ^ one can transform this last term into 
4G ^ 2\-)^ ^ 


X 


E 

p>0 


-FP 

^2n-p+l 




(2n-p+Z) 


|y|<^ 


T (i)-M(r*^)(y,t- |y|/c). (5.20) 


The superscript {p + l) on the yrfunction refers to the ^-differentiation. It is straightforward 
to show, using the fact that 7 /( 2 :) = (—-|- 1)!!P;(2:) is directly related to the Legendre 
polynomial, that 


(p^9(i) ^ ^ (2/-H)!!(/-Hp)! 


(5.21) 


2Pp\{l — p)\ 

Since Al(r*^) is singular at the origin (but regular at inhnity), and because of the explicit 
factor B brought about by the derivative 9i(|y|'®), the integral in Eq. ( ^.201) ranges over 
|y| < TZ (and even |y| < e where e is an arbitrary small number). We can then expand 
7 Vl(r*^)(y,t — |y|/c) when c —>■ -l-cxo. Furthermore we can change the integration over 
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|y| < TZ into an integration over |y| > TZ hj simply changing the sign in front of the 
integral. Indeed, this comes from a technical lemma, which plays an important role in Refs. 
|, |i0[| ; see before Eq. ( |C6D in Appendix y, and the proof given in Ref. [^]. Thus, 


AG ^ 2^{-y , 

■ c 4 ( 2 / + l)!^ 

n,l,k>0 ^ ' 


■FP 


X V_^ 

^ c2"-!>+'+‘S;0 




(2n—p+/+/c) 


c?V5.(|yrK^|yr-^-' T (i)A^(0(y,t) • (5.22) 

k\ 


By changing the label k into 2k + p — I and 2k + 1 + p — I, in order to cover odd and even 
numbers, we are able to write the previous expression in terms of some sums of real numbers 
indexed by p, i.e. 


AG ^ 2‘(-y 


n,/,/c>0 


(p-l + 2k)! 


xFP 

B=0 


(2n-l-2k) 

o!Va.(|y|®)*L|y|“-'-'A<(T"')(y,i) 


|y|>^ 


AG ^ 2^(-y (-] 


' (-ir'T(i) 


n,/,/c>0 


xFP 

B=0 


p=0 

iB 


l + 2k + l)\ 


(2n+2fc+l) 


d^y mnnLlyr-^ Mir^ (y,t). 


(5.23) 


|y|>'?^ 


The sums in curly brackets are found to be explicit expressions depending on k and I and 
some factorial combinations : 


Vk>t + 1, 


( 2 / + 1 )!! 


p'^o (P-^ + 2^)! 2>^k\{2k -21- 1)!! ’ 


v,>, ±^rG:\i) 


( 2 / + 1 )!! 


(p - / + 2A; + 1)! 2^-\k - l)\{2k + 1)!! ’ 

^k<i, A (-ir'T(i) j + 2*^-1)!! 


^ (p-/ + 2fc)! 

v,<,_, ^ (-ir-T(i) 

^ (p - / + 2fc + 1)! 

p=0 ' ■' 


2m 


(5.24) 

(5.25) 

(5.26) 

(5.27) 


Thanks to these formulas, one can transform Eq. (|5.23|) into 
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(2n+2fc) 


J\y\>n 
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(2n+2fc) 
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^V5*(|yr)nL|y|“-M(n(y,t) 


A-«(fL) 
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xFP 
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'|y|>^ 


d^ydmnnL\yr-^M{T^n (y,t). 


In the latter expression we can recognize 
Vfc > /, A-*^(aL(|y|~^)) = 
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{2k - 21 - l)\\2m 


nL\y 


2k-l-l 
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2k-l-l 


A ^(^l) = 


2^fc! 
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SO that we obtain 
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(2n+2l+l) 


(5.28) 
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E /n~,\ t, ^~"fe)FP / 4VSi(|yr)9L / d27i(z)A<(T“')(y.«-5|y|/c) 
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n,/,fc>0 






__ (2n+2/c) 

ASi(|y|pA-‘(4(|yr‘))X(i^’'))(y.«) 


A E(z4A-.(i.)FP 


C’ — (2/ + 1)! 

n,l,k>0 ^ ’ 


B=0 


{2n+2k+2l+l) 

d^yd,{\yn^-'^iyL)M{rn{y,t) . (5.32) 


|y|>'?^ 


After replacing Eq. (|5.32|) in Eq. (|5.19|) , at long last we hnd 


= 0 . (5.33) 

In this way, we have checked that the post-Newtonian metric, found by matching as a 
dehnite functional of the stress-energy pseudo-tensor satishes the harmonic-coordinate 
condition as a consequence of the conservation of this pseudo-tensor. 
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APPENDIX A: NEAR-ZONE EXPANSION OF THE RETARDED INTEGRAL 


This appendix, provided here for completeness, is an extended, and also somewhat sim- 
plihed, version of the derivation given in Appendix A of Ref. |^. We are interested in 
source functions, say Af(r)(x, t), having the form of an exterior multipole-moment decom¬ 
position, valid outside the compact-support domain of the source. We employ the same 
notation as in Section |II A| (except that we do not write the space-time indices) : r denotes 
the pseudo-tensor of the source; notably we have M.{t) = where A is the grav¬ 

itational source term. The two basic properties of the function Al(r)(x, f) are that it is 
smooth on deprived from the spatial origin r = 0 : 


Wl(r)(x,f)eC''-(M3x 


(Al) 


and that it admits a near-zone expansion, when r —0 (with t = const), having the appro¬ 
priate structure [c/Eq. (|2.14|) ] : i.e., VA G N, 


M (r) (x, t) = hLr“(ln rYGL,a,p(,'t) + , 


(A2) 


where a G Z with a < N and p G N. 
formal (inhnite) near-zone expansion. 


Like in Section IllH we denote with an overline the 


Wl(r)(x,f)= '^hLr‘"{\nryGL,a,p{t) ■ (A3) 

It is very important to make the distinction between Ai (r) and its formal near-zone expan¬ 
sion Here we shall investigate the retarded integral of the product (r)(x, t), 

where R G C, by means of analytic continuation (we pose rg = 1 in this Appendix). For this 
task we assume at hrst that the real part of B is large enough so as to “kill” the divergencies, 
when r —>■ 0, of the expansion ( |A2|) , so that the retarded integral is initially well-defined. 
Therefore, rigorously speaking, we are allowed to do this only if there exists a hnite maximal 
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divergency, i.e. some Omin < a in Eq. (^) with finite amin £ We have seen in Section 
Eg that such maximal divergency exists at any given post-Minkowskian order m, but no 
longer exists for the full post-Minkowskian series because amin(’^) —oo when m —>■ -|-cxo. 
The consequence is that the analytic continuation is in principle justified only at a given 
finite post-Minkowskian order. But, as explained in Section |11 A| , we sum up systematically 
all the post-Minkowskian results. In this way we are entitled to proceed as we do below; 
simply we have to remember that the end result will be a priori true only in a sense of 
formal post-Minkowskian expansions. 

We decompose the source term into multipoles according to 


+ 00 


Af(r)(x,t) = '^fiLaLir.t) , 


(A4) 


1=0 


where the cr/,’s are STF functions in L = ■ ■ - b. The inverse formula is 


= - - - / —nLM(T)(x,t) , 


(AS) 


where dQ is the solid-angle element around the unit vector rii = rd = x^jr. Then the 
expression of the retarded integral, in a sense of analytic continuation in 5, is given by the 
following explicit formula, obtained in Ref. [see Eqs. (6.3)-(6.5) there] : 


|•t—r \ jdB ( t—r—s \ pB ( t+r—s \ 

□Re\ [r^M (r )ds d, ^ ^ 

1=0 ^ 

(we pose c = 1 and ro = 1 in this Appendix), where the function R^^p, s) reads 


(A6) 


Rf(p, s)=p*y dx ^ x^aL{x,x + s). 


(A7) 


l\ \x^ 

Following the same procedure as in Eqs. (A6)-(A7) in Ref. [^ we are allowed to re-write 
the expression (|A6|) into the alternative form 
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°Re\ [r^-M(r)(x,t)] 


~|“00 n'f' 

dudL‘{ ^Rl 

1=0 d-r 


1 / M + r 


r 


,t — u 


+ O0 


--y 

A'tt 


47r l\ 
1=0 


-Bl -nf{t + r) 


2r 


(A8) 


The “anti-symmetric” wave is parametrized by 7?.f(t) which is related to the function 
Rf{p:s) by 


nf(t) = S7r(-y+H\ dsRf(^ 


t — s 


(A9) 


Inserting Eq. o, and performing some change of variables, we obtain 


Kit) = 


Anil ^ / 

^ J dxx ++ y dz-fi{z)aLixR- zx) , 


(AlO) 


and, using the relation (^), and considering the variable x as the norm of x G M^, we 
further get 


TZf{t)= I d^x / dz'yi{z) M{T){-K,t — z\x.\) . 


iB 


^+oo 


(All) 


In these expressions the function jilz) is defined by 


7i(.-) = (-)“2E^(z"-l)'. (A12) 

where the particular /-dependent factor has been chosen in such a way that the integral is 
normalized to one in the following sense (see Ref. ). Considering first that I is a complex 
number such that —1 < 3?(/) < —1/2 we can compute the integral of 'yi{z) by means of the 
Euler T-function, with the result 


f + OO 


dz-fi{z) = 2(-) 


,+, r(2/ + 2)r(-2/-i) 
r(/ + i)r(-/) 


(A13) 
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The right-hand-side of this equation can be analytically continued to all values I G C except 
half-integer values, and is found to be equal to one when I is an integer : 


+ 00 

dz-fi{z) = l (/gN). (A14) 

Next, let us treat the hrst term in the right-hand-side of Eq. (|A8|) , say 

= . (A15) 

This term is a particular solution of the d’Alembertian equation = r^Ai{T) [since the 
second term in Eq. (A8|) is a source-free solution]. We shall prove that the (formal) near¬ 
zone expansion of that term, i.e. J^(x, t), is given by the integral of the “instantaneous” 
potentials acting on the near-zone expansion of the source term, i.e. r^Ad(r)(x, t). For any 
of the terms composing the multipolar source r^Af(r) [see Eq. (|A3|)], we hrst dehne 



A ^ [hLr^+“(lnr)PGL,a,p(t)] 


(dy 

hLr^+“+2GL,a,p(t) 

\dB) 

_{B + a + 2 - 1){B + a + 3 + 1)_ 


(A16) 


(this being justihed by the fact that one gets an identity by applying A on both sides). 
Clearly the previous formula can be iterated and so we can dehne the operator A“^“^ = 
(A-i)fc+i, applied on each separate terms in Eq. (|A3|) and therefore on the complete series 
(r)(x, t). From this we obtain the instantaneous-potentials operator, as the formal 
expansion series 




r^Al(r)(x, f) 



r^Af(r)(x, f) 


(A17) 


Notice that this operator contains only some even powers of 1/c. An important point 
for our purpose is that [r'®Af (r)] is proportional to the regularization factor r'®; and 
it evidently satishes □(X“^ [r®Af (r)J) = r'®Al(r). On the other hand, we have also the 
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equation = r'®A^(r), which comes from applying the overline operation onto = 
r^M.{T). This shows that r^M(T) and must differ by a solution of the homogeneous 
equation, hence there should exist some functions C^(t) and D^{t) such that 


+ 00 




C^(t — r) + Df(t + r) 


(A18) 


z=o 


Note that the dependence on B of the second term is “hidden” inside the functions 
and D^. Let us now prove that in fact the latter functions must be zero. This is a simple 
consequence of the expression ( [A7| ) for the function s), from which we deduce that 

the expansion when p ^ 0 of this function is proportional to in fact, it has the structure 
Rl (p. s) ~ ^ pS+6(inp)<?, vi^hen p ^ 0. From this knowledge, we easily hnd that the near¬ 
zone expansion of is proportional to the factor . Since, as we have remarked, this is 


also the case of the hrst term in Eq. ( A18 ), X ^ [r'®Ad(r)j, and since it is impossible that 
(the near-zone expansion of) the second term in Eq. ( |A18|) be itself proportional to - the 


B’s affect only the functions and Df but not the structure of the near-zone expansion 
- we conclude that and are identically zero. Hence we have proved 


J^(x, t) = X ^ 


r^M(T)(x,t) 


(A19) 


It suffices now to apply the overline operation (i.e., to take the near-zone expansion) onto 
Eq. (|A8|) to get our hnal result. 


°Ret lr^M(T)(x,t)] =x 


-1 




-«?(*->■)-Kf(i + r) 

r^XT(r)(x, t) 

dTT^ n 

1=0 

2 r 


(A20) 

where we recall that the function (t) has been given by Eq. (|A11|) - [The formula used 


Section [II B| results from applying the hnite part operation FP .] 
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APPENDIX B: THE GENERALIZED POISSON OPERATOR 


In Appendix ^ we have been interested in source functions of the multipolar type 
AI(r)(x, f), which are smooth in x M and possess a near-zone expansion of the type 
(|A3|). In the present Appendix we consider some source functions of the post-Newtonian 
type r(x, f). These are supposed to be smooth all over 


r(x,t) e , 


and to admit a far-zone expansion with structure (VA G N) 


(Bl) 


r(x,f) = ^ (InrfGi,a,p(t) S'Ar(x,f) , 


(B2) 


where a G Z, with —N < a, and p G M. The remainder-term is S'Ar(x, t) = O (l/r^) when 
r —>■ -l-cxo with t = const. 

Let us consider some B E C, and a radius 7?. G M with 7^ > 0. We define two inte¬ 
grals, corresponding to a split of the Poisson integral between “near-zone” and “far-zone” 
contributions, separated by the radius TZ : 


= -S I 


d^y 
|x- y| 
d^y 


|y|^7"(y,^) • 


(B3) 

(B4) 


'\y\>Tl |X-y| 

The i7-dependent regularization factor is |y|'® = (|y|/ro)'®. It is easily checked that the near¬ 
zone integral /^(x, f) is well-dehned (convergent) when ^{B) > —3 and that the far-zone 
one /^(x, f) is well-dehned when 3fJ(i7) < —Omax — 2, where Omax is the maximal power of r 


in the expansion ([B^). So we have to assume at this stage the existence of some maximal 
divergency corresponding to some power Omax- Strictly speaking, our present investigation 
is thus valid only at some hnite post-Newtonian order. But, in fine, we sum up the results, 
and we consider the complete post-Newtonian series to hold true in a formal sense. 
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We want first to check that the integrals (|B3|) and (p^) can be analytically continned 
down to a neighbonrhood of B = 0 (except at the valne B = 0 itself), let say in the open 
domain dehned by 0 < \B\ < e (where e < 1). There is no problem with the near¬ 
zone integral J^(x, t) which is clearly convergent all over B^ and even at the valne B = 0. 
Concerning the far-zone integral />(x, t) we replace the fnnction r inside the integrand by 
its far-zone expansion (p2|) : 


= {^^n(y)|yr(ln|y|)^GL,a,p(t) + >gjv(y,t)| . (B5) 

When N is large enongh the contribntion dne to the remainder Sn is convergent all over B^ 
and at B = 0, with evidently the valne at B = 0 given by 


f d^y\y\^ 

l\y\>n |x-y| 


-SN{y,t) = 


d^y 


-SN{y,t) + 0{B) . 


(B6) 


\y\>Tl |X-y| 

Thns we need only to deal with the other contribntions, which consist of a hnite snm of 
terms, say 


Y, f ^^fti(y)lyr(in|yl)'’- (B7) 

J\y\>Tl |X y| 

Let ns snppose that the held point x lies inside the far-zone domain, i.e. TZ < |x|. We 
distinguish the two cases where |y| < |x| and |x| < |y|. For each of these two cases 
we substitute into the integrals the appropriate multipolar expansion of the factor 
for instance when |y| < |x|. This leads, after performing the 

integration over the angles, to some series of radial integrals having the structure (ignoring 
some unimportant factors) 

/ d\y\\y\^+<^+^+\\n\y\r + / ^lyl |y|)" • (B8) 

|x| Jn j|x| 
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When \x.\ < TZ the reasoning is the same but simply one ignores the first term in Eq. ( P8|) 
and take TZ as lower bound in the second term. Computing each of these integrals we find 


E 


XL 

( dV 

X 

X 2*+l 

KdBj 





E 


XL 


_d_ 

dB 


p r 


\B+a-l+2 


B + a-l + 2 


(B9) 


B+a+l+3 

Each of these terms clearly admits an analytic continuation for any B E and in fact 
for any B E C except at integer values. Furthermore we see from that expression that the 
function will admit a Laurent expansion when B ^ 0, with in general some multiple poles 
[coming from the differentiation {d/dBy of simple poles ~ 1/B]. Hence our statement. 

It is clear that the Laplacians of the two integrals and /> satisfy, in the domains of 
the complex plane where these functions were initially valid : 


3?(H) > -3 
3?(H) < -a^ax - 2 


A/^(x, t) = Y{TZ— |x|) |x|'® r(x, t) , 
A/^(x, t) =Y (|x| — TZ) |x|'® r(x, t) , 


(BIO) 

(Bll) 


where Y denotes the Heaviside step-function. Therefore, if we define for any B E Be the 
object 


/^(x, t) = (x, t) + analytic continuation < /> (x, t) \ , (B12) 

we hnd that it necessarily satisfies, for any B E Be, the H-dependent Poisson equation 

AJ'®(x, f) = Ixl"® r(x, f) . (B13) 

On the other hand, we have learned from Eq. ( |B9| ) that admits when H —> 0 a Laurent 
expansion involving (in general) simple and multiple poles. Now the key idea, as we shall 
prove, is that the finite part, or coefficient of the zero-th power of B in the latter Laurent 
expansion, represents a particular solution of the Poisson equation that we want to solve. 
Let the Laurent expansion of be 
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(B14) 


H-oo 

/^(x,t)= ^ 4(x,t)5^, 

h—k • 

— '^min 

where fcmin G and where the coefficients ik depend on the held point (x, t). By applying 
the Laplacian operator onto both sides of Eq. (|B14|) , and using the result ([B13|) together 
with the Taylor expansion of the regularization factor |x|^, we arrive at 

/Cmin <k < -I Aik = 0 , (B15) 

, ^ (Inlxl)^ , , 

k>0 A4 = r . B16 

k\ 

Thus, the case fc = 0 shows that the hnite-part coefficient in the expansion (|B14|) , namely 
io, is a particular solution of the required equation : Azq = T. We shall now forget about 
the intermediate name zq, and denote, from now on, the latter solution by A~^t = zq, or, in 
more explicit terms. 


A“ir(x, f) = FP A ^ Ixl^rfx, t) 

^ ^ B=0 L , 


(B17) 


where A“^ refers to the standard Poisson integral, and the hnite-part symbol FP means the 
previous operations of considering the Laurent expansion when B ^ 0, and picking up the 
hnite-part coefficient. Thus, we have proved that A[A“^r] = r, so the generalized inverse 
Poisson operator A“^ dehnes a particular solution of the Poisson equation, which has, by 
construction, none of the problems of divergencies of Poisson integrals which have so much 
plagued the standard post-Newtonian approximation |^-[T^ . 


Finally let us prove that our generalized solution A'^r owns the same properties (pTI)- 
( |B2|) as the corresponding source r. This verihcation is important because it will allow us to 
iterate any number of times the operator A“^, and to obtain the post-Newtonian expansion 
up to any post-Newtonian order. The main problem amounts to prove that A“^r admits 
the same type of expansion at inhnity |x| —> +cxo as in Eq. (P2|). To do this we consider 
again the same split into near-zone and far-zone contributions ; A“^r = /<-!-/>, where 
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/<(x,t) = - —FP / 
4vrs=o 

/>(x,t) = -^FP 


d?y 

|x-y| 

d?y 


|ylMy,^) ■ 


(B18) 

(B19) 


47rs=o |x - y 

The near-zone integral admits an expansion at infinity which is of the required type. Indeed, 
because the integrand is of compact support, |y| < 7^, we can replace in it the factor 
by its expansion ^ ^^*4 integrate term by term. So we have, 'iN G N, 


n 

N-l 




dvr /! 
1=0 


-dr 


- ) FP 

r / B=o 


|y|<^ 


d^y\yry^T{y,t) + o 


„N 


(B20) 


The right-hand-side has indeed the same structure as in Eq. (|B2|) . The treatment of the 
far-zone integral is more delicate. We proceed in a way similar to what was done in Eqs. 
( ^SD - fpOf) . Namely we replace into it the source r by its expansion given by Eq. (|B2D . This 
yields (the finite part of) Eq. (pSf) , that for convenience we reproduce here : 


4(x,t) = -^FP / ^y|yl |^^^(y)|y|°(in|y|)PG'^ (t)+ gjv(y,t)| . (B21) 

47rB=o |x - y| J 

There is a contribution of the remainder and a finite sum of terms with known structure. 
The remainder contribution is simply given by the value at B = 0 which has been written 
in the right-hand-side of Eq. (|B(j|) . Let us write this term in the form 


'\y\>n |x-y 


Af-4 


■SN{y,t) = 




1=0 


d^yy^Sniy, t) + ^w-2(x, t) , (B22) 


7|y|>7^ 


where we introduced the — 4 first terms of the multipolar expansion of when r = 


-l-cxo, and where 


rAr_2(x,f) = 


'|y|>^ 


d^y 


N-i 


-E 

/=o 


/! 


-y^d. 


SN{y,t) . 


(B23) 


The maximal order iV — 4 of the expansion is chosen in such a way that all the terms in 
Eq. (|B23|) are given by convergent integrals at infinity, owing to the fact that the remainder 
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satisfies Sn = 0{l/r^). Now we prove that T 7 V- 2 , defined by Eq. (P23|) , is also a remainder 
in the sense that T;v _2 = 0{lnr We split Ti ^_2 into two integrals, a near-zone 
integral corresponding to the integration range |y| G ]7^, |x| [, and a far-zone one 
corresponding to |y| G ] |x|, -|-cxd[. In the near-zone integral we can use the bound 

< • (B24) 

where Cn is a constant. On the other hand, because Sn = 0{l/r^), there is also a constant 
An, depending on the value of 71, such that the following majoration holds : 



|57v(y,f)| < 



Replacing these results into the near-zone integral we get 


(B25) 


() [ < te 111 j , (b26) 

In the far-zone integral, we can no longer apply the bound ( |B24|) but still we can employ 
the majoration ([B25|). Then we can easily show the inequality (in which |y| = |x|A) 


T7V-2(x,i) < dTT 


A 


N 


r’+oo 


\N-2 


dX 


N-4 


(2/ - 1)- y 


z=o 


n 


(B27) 


The integral is convergent. At last, from Eqs. (|B2b| ) and (B27) we have proved that 
Tn- 2 = 0{\nr/r^ ^). Still it remains to show that the hnite sum of terms in Eq. ([B21|), i.e. 


besides the remainder, admits some expansions of the required structure. But this follows 
from applying the hnite part operation FP onto the result (p^) , which tells us immediately 

B=0 

that we have an expansion of the correct type ~ hi(x)|x|“(ln Ixl)”?. 
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APPENDIX C: FAR ZONE EXPANSION OF THE POISSON INTEGRAL 


Thanks to the investigation in Appendix the far-zone (or mnltipolar) expansion of the 
object A~^[r] happens to be workable. Recall that controlling the far-zone expansion of the 
post-Newtonian held is fnndamental since it is at the basis of the matching. The operation 
of taking the far-zone expansion is denoted A4 when applied on post-Newtonian objects (see 
Section [III B|) . We therefore want to determine the expression of Ad(A“^[r]). That is, we 
want to relate it to the expansion of the corresponding sonrce, which has the same strnctnre 
as in Eq. (|3.22|) : 


M(T)(x,t)= Y^nLr‘^(lnr)PGL,a,p(i) ■ (Cl) 

By the matching eqnation we know that this far-zone expansion is identical with the near¬ 
zone expansion of the external held [see e.g. Eq. (|A3|) ]. Let ns hrst apply M onto A“^[r] 
as expressed as a snm of near-zone and far-zone contribntions, 


V((A->lT])=Vf(/<)+>((/>), (C2) 

where /< and /> are dehned by Eqs. (|B18| )-( pT^ . The near-zone integral is qnite easy 
to work with. Indeed, from Eq. ([B20|) we see that its expansion when r = |x| —>■ -|-cx) 
is obtained by expanding the factor l/|x — y| inside the integrand. Therefore, the inhnite 
far-zone expansion (withont remainder) reads 

|y|-A. (^) T(y.t) . (C3) 

in which we denote 


M 


|x-y| 


+ 00 

E 

1=0 


-VlOl ^ 


(C4) 
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On the other hand, the far-zone expansion of the far-zone integral /> has been obtained 
in Eq. ( |B21| )- (|E^ ), where we found that it comes from replacing the source term by its 
far-zone expansion [indeed, when TZ is large enough, the integration ranges over the domain 
of validity of the far-zone expansion]. So the inhnite far-zone expansion of that term is given 
by 


M{I,) =-^FP f p^M(ny,t)), 

^ 47rB=o |x-y| V V ’ 


(C5) 


where the integrand contains the expansion of the source given by Eq. m- Now let us use 
a technical lemma which is quite important in the present formalism, and has already played 
a crucial role in Refs. P^,|40]. This lemma is based on the remark that any radial integral 


of the type fg °° (i|y| |y|'®+“(ln |y|)P, where R G C and a and p are arbitrary real numbers, 
is identically zero by analytic continuation in B. See Ref. for the proof. Our useful 
lemma, that is trivial to relate to the previous remark (after performing the integration over 
angles), is 


FP I dV \y\BM M (r(y, f)) = 0 . (C6) 

The point here is that the integral ranges over the complete three-dimensional space 
Now we have the “numerical” equalities Af when |y| < |x| and A1 (r) = r 

when |y| > a, where a is the radius of the compact-support source. From this we deduce 
that as soon as TZ > a, what we can always assume right from the beginning, and |x| > TZ, 
which is not a problem because we are considering the limit |x| —>■ -|-cxo, we have the identity 

FP f xfr(y,i)') + FP f <i“y |y|®M ( 7 — !— 7 )’'(y.t) = 0 . (C7) 

By means of that identity we can obtain the requested form of the far-zone expansion as 


+M f , _ | j 


• (C8) 
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In this particular form we see that the Al-operator is distributed on the two terms like a 
derivative operator would be. In the hrst term we recognize the action of the generalized 
Poisson integral. Actually this Poisson operator has been dehned in Appendix ^ when 
acting on a near-zone expansion of the type (|A3|) , but by matching that expansion is the 
same as the present far-zone expansion, so the dehnition is rigorously the same. Finally we 
can re-write Eq. ([C^) into the alternative form 


^ ^ __ -| +00 f \l /» 

M \m{t)] - t-E fP / i'y ly|'’»iT(y, t), 

" 1=0 ' 


(C9) 


which constitutes the main result of this Appendix. Notice that Eq. ( |C9|) is in agreement 
with the multipole expansion of the retarded integral as given by Eq. (3.11)-(3.12) in Ref. 
ig , when specialized to the static case where there is no dependence on time. 

Next we derive the analogous result concerning the operator of the “instantaneous” 


potentials 



k=0 


(CIO) 


We iterate k + 1 times the result (|Cg). There is no problem for doing this; the only point 
is that we use in a repeated way the easily checked formula telling that we are allowed to 
“operate by parts” the Poisson integral in the way 


FP j d^z\ z\^zlA~^t= FF j d^y\y\^A-^[yL]T 


This formula is a consequence of the fact that FP j d^z \z 


B zl — 


3 = 0 '' ' ' i^-yi 


(Cll) 
-47rA-i [Dl] = 


-27r |,a2 


21+3 


yyyL', see Eq. (4.10) in Ref. [^. Therefore we arrive at 


A4 [r]j = 


(C12) 


1 


+00 / \/ ^ 


Z=0 i=0 


FP d^y\y\^/\^-^[yL]T{y,t) , 
£>=0 / 
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and from this it is very simple to derive the requested expression concerning We obtain 


=X-i[Af(r)] 


(C13) 


+ CXD 


\l +°° 


+ 00 


1 


Z=0 i=0 k=i 


This expression, though completely explicit, does not constitute our hnal form. Because the 
“instantaneous” solution is a particular solution of the d’Alembertian equation, it must be 


possible to re-express the second term in Eq. (|C13|) as a combination of some source-free 
retarded and advanced multipolar waves. To see this we notice that 


A-* 


dhir 


= dr 


„2i-l 


( 2 !) 


(C14) 


which shows that the latter homogeneous solution is actually one of the symmetric type, i.e. 
retarded plus advanced. Namely we can re-write Eq. (|C13|) into the form 


I ( \i 




dvr /! 
1=0 


(-) ^^ J ^Lit - r/c) + TL{t -t- r/c) 


2r 


(CIS) 


where the overline notation means taking the Taylor expansion of the symmetric wave when 
the retardation r/c —> 0 [the result is displayed in Eq. (|3.24|) ]. Actually, this overline 


notation is somewhat misleading, because, in keeping with the real meaning of the result 
( CIS ), one should a posteriori interpret the latter Taylor expansion as a far-zone (singular) 
expansion when r —> -|-cxo. However, in view of the matching, it is more fruitful to employ 
the same overline notation as for the expansion of the anti-symmetric waves occuring in 
the near-zone metric - indeed the matching is simply interested at identifying together some 
asymptotic expansions which are of the same form. The “multipole-moment” function tFiit) 
in Eq. (|C15|) is given by 


+00 .. „ 

^L{t) = t) 

j=o 


(CI6) 
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Finally let us find an alternative, more compact, form for this result. We introduce the 
/-dependent function 


^ {21 + 1)!! _ 2y 

2i+ii\ ^ ’ 


(C17) 


whose integral is normalized to one : dz6i{z) = 1. One can readily show that 


[vl] = |y, (CIS) 

which permits to express the function in a form where the post-Newtonian series is 
formally re-summed as 


^l(<) = FP j <i’x|y|V r 


dz6i{z)r{y,t± z\y\/c) . 


(C19) 


Under this form we recognize the multipole-moment function introduced in Eq. (3.14) in 
Ref. (the function remains unchanged by taking either sign ± in the time argument of 
r). This result permits to fully determine the exterior multipolar held by matching, and 
to recover the expression already obtained in Ref. by means of a somewhat different 
method. 


57 







REFERENCES 


[1] A. Einstein, Sitzber. Preuss. Akad. Wiss. (Berlin), 688 (1916). 

[2] J.Droste,Versl. K. Akad. Wet. Amsterdam 25, 460 (1916). 

[3] W. DeSitter, Mon. Not. R. A. S. 76, 699 and 77, 155 (1916). 

[4] A. Einstein, L. Infeld and B. Hoffmann, Ann. Math. 39, 65 (1938). 

[5] H.A. Lorentz and J.Droste, Versl. K. Akad. Wet. Amsterdam 26, 392 and 649 (1917); 
in the collected papers of H.A. Lorentz, vol. 5, The Hague, Nijhoff (1937). 

[6] V. Fock, J. Phys. (U.S.S.R.) 1, 81 (1939). 

[7] V.A. Fock, Theory of Space, Time and Gravitation, Pergamon, London (1959). 

[8] S. Chandrasekhar, Astrophys. J. 142, 1488 (1965). 

[9] S. Chandrasekhar and Y. Nutku, Astrophys. J. 158, 55 (1969). 

[10] S. Chandrasekhar and F.P. Esposito, Astrophys. J. 160, 153 (1970). 

[11] J.L. Anderson and T.C. DeCanio, Gen. Relat. Grav. 6, 197 (1975). 

[12] J. Ehlers, in the Proc. of the International school of Relativistic Astrophysics, J. Ehlers 
(ed.), Erice, MPLMiinchen (1977). 

[13] J. Ehlers, Ann. N.Y. Acad. Sci. 336, 279 (1980). 

[14] G.D. Kerlick, Gen. Rel. Grav. 12, 467 (1980). 

[15] G.D. Kerlick, Gen. Rel. Grav. 12, 521 (1980). 

[16] A. Caporali, Nuovo Cimento 61B, 181 (1981). 

[17] A. Papapetrou and B. Linet, Gen. Relat. Grav. 13, 335 (1981). 

[18] R. Brener and E. Rudolph, Gen. Rel. Grav. 13, 777 (1981). 


58 



[19] R. Breuer and E. Rudolph, Gen. Rel. Grav. 14, 181 (1982). 

[20] T. Damour and N. Deruelle, Phys. Lett. 87A, 81 (1981). 

[21] T. Damour, G. R. Acad. Sc. Paris 294, 1355 (1982). 

[22] T. Damour, in Gravitational Radiation, N. Deruelle and T. Piran (eds.), North-Holland 
Gompany (1983). 

[23] L. Blanche!, G. Faye and B. Ponsot, Phys. Rev. D58, 124002 (1998). 

[24] P. Jaranowski and G. Schafer, Phys. Rev. D57, 7274 (1998). 

[25] T. Damour, P. Jaranowski and G. Schafer, Phys. Rev. D62, 021501R (2000); Erratum 
Phys. Rev. D63, 029903 (2001). 

[26] L. Blanche! and G. Faye, Phys. Lett. A 271, 58 (2000). 

[27] L. Blanche! and G. Faye, Phys. Rev. D63, 062005 (2001). 

[28] V. G. de Andrade, L. Blanche! and G. Faye, Glass. Quantum Grav. 18, 753 (2001). 

[29] W.L. Burke and K.S. Thorne, in Relativity, M. Garmeli et al. (eds). Plenum Press, New 
York, p. 208-209 (1970). 

[30] W.L. Burke, J. Math. Phys. 12, 401 (1971). 

[31] L. Blanche! and T. Damour, Phys. Rev. D37, 1411 (1988). 

[32] L. Blanche!, Phys. Rev. D47, 4392 (1993). 

[33] L. Blanche!, Phys. Rev. D55, 714 (1997). 

[34] B.R. Iyer and G.M. Will, Phys. Rev. Lett. 70, 113 (1993). 

[35] B.R. Iyer and G.M. Will, Phys. Rev. D52, 6882 (1995). 

[36] L. Blanche! and T. Damour, Philos. Trans. R. Soc. London A320, 379 (1986). 


59 



[37] L. Blanchet and T. Damour, Ann. Inst. Henri Poincare A 50, 377 (1989). 

[38] T. Damour and B.R. Iyer, Ann. Inst. Henri Poincare A 54, 115 (1991). 

[39] L. Blanchet, Phys. Rev. D 51, 2559 (1995). 

[40] L. Blanchet, Class. Quantum. Grav. 15, 1971 (1998). 

[41] L. Blanchet, T. Damour and B.R. Iyer, Phys. Rev. D 51, 5360 (1995). 

[42] C.M. Will and A.G. Wiseman, Phys. Rev. D 54, 4813 (1996). 

[43] L. Blanchet, Phys. Rev. D 54, 1417 (1996). 

[44] L. Blanchet, Class. Quantum. Grav. 15, 113 (1998). 

[45] L. Blanchet, B.R. Iyer and B. Joguet, Phys. Rev. D 65, 064005 (2002). 

[46] A.D. Rendall, Proc. R. Soc. London A438, 341 (1992). 

[47] T. Futamase and B. F. Schutz, Phys. Rev. D 28, 2363 (1983). 

[48] T. Futamase, Phys. Rev. D 28, 2373 (1983). 


[49] Our notation for STF tensors is the following. L = iii 2 ■ ■ - ii denotes a multi-index, made 
of I (spatial) indices. When summing over multi-indices we never write the I summations 
over the I indices A, • ■ ■ ,ii ranging from 1 to 3. The STF product of unit vectors 
rii = rf = xQr is denoted hL=STF(nL), where ul is a short-hand for ■ ■ -riii. For 
instance hij = Uirij — ^Sij. Similarly we denote xl = Xi^- ■ -x^ = r^rii and fL=STF(a:L). 
The derivative operator Ol is a short-hand for di^ - ■ ■ di^, and we have dL=STF{dL). For 
instance dij = dij — A. More generally, a function Fl is said to be STF with respect 
to the I indices composing L if and only if, for any pair of indices ip, ig G L, we have 
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comes from the non-stationary (or radiative) part of the held, which according to our 
assumption of stationarity in the past is zero when t < —T, and for which TZ^ [t) is 
perfectly well-dehned. For simplicity in the notation we do not indicate that TZ^ {t) 
should be computed only from the “radiative” part of the source term 

[51] We want to prove that the radial integral /o^°°d|y| |y|^’''“(ln |y|)^ is zero by analytic 
continuation (Vi? G C). First we can get rid of the logarithms by considering some 
repeated differentiations with respect to B; thus we need only to consider the simpler 
integral /g^°°d|y| |y|'®’'"“. We split the integral into a near-zone integral d\y\ |y|'®’''“ 
and a far-zone one d\y\ |y|^’'"“, where IZ is some constant radius. When 3?(i?) is a 
large enough positive number, the value of the near-zone integral is 7^^’''“’''^/(i?-|-a-|-l), 
while when 3fJ(i?) is a large negative number, the far-zone integral reads the opposite, 
—7Z^^°‘^^/{B + a + l). Both obtained values represent the unique analytic continuations 
of the near-zone and far-zone integrals for any i? G C except —a — 1. The complete 
integral /Q^°°(i|y| lyl^"*"^ is equal to the sum of these analytic continuations, and is 
therefore identically zero (Vi? G C, including the value —a — 1). 
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